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The  purpose  of  this  study  is  to  investigate  and,  if  possible,  establish  the 
feasibility  of  producing  artificial. geomagnet  ic  micropulsations  by  using  a  large 
horizontal  current  loop  on  the  ground.  Most  natural  micropulsations  appear  to  be 
produced  by  hydromagnetic  (hm)  waves  propagating  in  the  ionosphere  and/or  the 
magnetosphere,  and  the  term  "artificial  micropulsations''  is  used  to  mean  artifi¬ 
cially  produced  hm  waves  with  propagation . characteristics  similar  to  those  of  the 
natural  signals. 


The  critical  quantity  in  the  generation  process  for  the.hm  waves  is  the  amplitude 
of  the  magnetic  :ield  variation  that  can  be  produced  in  the  ionospheric  E-region 
by  the  ground-based  loop.  Currents  are  induced  in  the  conducting  earth  beneath 
the  loop  whenever  its  current  is  varied,  and  these  currents  have  a  magnetic  field 
opposing  the  magnetic  field  of  the  loop.  If  either  the  conductivity  of  the  earth 
or  the  variation  frequency  is  too  high  the  field  at  E-region  heights  may  be 
reduced  to  a  level  where  negligible  hm  wave  amplitudes  are  produced.  The  effect 
of  these  earth  currents  is  calculated  for  an  assumed  single- layer  earth  and  for 
a  wide  range  of  frequencies  (ltT2  -  102  Hz)  and  earth  conductivities  ,  (1<r  5~  4  mho/m). 
It  is  found  that  at  1  Hz  and  for  a  conductivity  of  10”^  mho/m  (representing  a 
typical  conductivity  for  the  earth's  crust.)  the  magnetic  field  of  the  loop  is 
reduced  by  about  18  percent.  At  lower  frequencies  the  reduction  in  magnetic  field 
is  sven  less,  Thus,  provided  the  loop  is  suitably  located,  the  effect  of  earth 

currents  should  not  be  very  great .  , 

__  _ _  (continued  on  page  2) 
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Expressions  for  the  two  hm  wave  modes  excited  in 
the  ionosphere  by  a  current  loop  on  the  ground  are 
derived  and  the  feasibility  of  constructing  large 
current  loops  is  considered.  It  is  concluded  that 
current  loops  with  magnetic  moments  as  large  as 
7  x  1013  amp>m2  are  practicable  and  that  the  cur¬ 
rent  in  these  loops  can  be  varied  at  micropulsa¬ 
tion  frequencies.  All  things  considered,  the 
total  magnetic  field  produced  in  the  JJ-region 
could  have  an  amplitude  as  large  as  0.  S  gamma, 
which  is  greater  than  most  natural  micropulsation 
amplitudes  observed  beneath  the  ionosphere. 

Preliminary  calculations  with  the  expressions 
for  the  two  hm  wave  modes  indicate  that  signifi¬ 
cant  hm  wave  amplitudes  can  be  produced  in  the 
ionosphere  above  the  larger  current  loops  con¬ 
sidered  in  this  study.  This  conclusion  is  in¬ 
dependently  confirmed  by  Greifinger  (1972).  Thus, 
we  conclude  that  the  artificial  generation  of 
ihicropulsations  is  feasible  and  that  a  large  cur¬ 
rent  loop  should  be  capable  oi  producing  observabli 
Pc  1  micropulsations  at  very  considerable  dis¬ 
tances  (greater  than  5000  km)  from  the  source.  '* 
Much  further  study  is  required  to  establish  all 
the  possibilities  of  the  artificial  generation 
process,  particularly  at  low  mi^ropulsation  fre¬ 
quencies  and  for  field-line  propagation.  However, 
a  number  of  important  uses  for  artificially- 
generated  micropulsations  are  immediately  apparent 
and  are  discussed  briefly. 
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ABSTRACT 


The  purpose  of  this  study  is  to  investigate  and,  if  possible, 
establish  the  feasibility  of  producing  artificial  geomagnetic  micro¬ 
pulsations  by  using  a  large  horizontal  current  loop  on  the  ground.  Most 
natural  micropulsations  appear  to  be  produced  by  hydromagnetic  (hm)  waves 
propagating  in  the  ionosphere  and/or  the  magnetosphere,  and  the  term 
"artificial  micropulsations"  is  used  to  mean  artificially- produced  hm 
waves  with  propagation  characteristics  similar  to  those  of  the  natural 
signals . 

The  critical  quantity  in  the  generation  process  for  the  hm  waves 
is  the  amplitude  of  the  magnetic  field  variation  that  can  be  produced 
in  the  ionospheric  E-region  by  the  ground-based  loop.  Currents  are 
induced  in  the  conducting  earth  beneath  the  loop  whenever  its  current 
is  varied,  and  these  currents  have  a  magnetic  field  opposing  the  magnetic 
field  of  the  loop.  If  either  the  conductivity  of  the  earth  or  the  variation 
frequency  is  too  high  the  field  at  E-region  heights  may  be  reduced  to  a  level 
where  negligible  hm  wave  amplitudes  are  produced.  The  effect  of  these 

earth  currents  is  calculated  for  an  assumed  single-layer  earth  and  for 

-2  2  • 
a  wide  range  of  frequencies  (10  -  10  Hz)  and  earth  conductivities 

-5 

(10  -  4  mho/m)  .  It.  is  found  that  at  1  Hz  and  for  a  conductivity  of 

10  4  mho/m  (representing  a  typical  conductivity  for  the  earth's  crust) 
the  magnetic  field  of  the  loop  is  reduced  by  about  18  percent.  At  lower 
frequencies  the  reduction  in  magnetic  field  is  even  less.  Thus,  provided 
the  loop  is  suitably  located,  the  effect  of  earth  currents  should  not  be 
very  great . 
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Expressions  for  the  two  hm  wave  modes  excited  in  the  ionosphere 

by  a  current  loop  on  the  ground  are  derived  and  the  feasibility  of 

constructing  large  current  loops  is  considered.  It  is  concluded  that 

13  2 

current  loops  with  magnetic  moments  as  large  as  7x10  amp.m  are 
practicable  and  that  the  current  in  these  loops  can  be  varied  at  micro¬ 
pulsation  frequencies.  All  things  considered ,  the  total  magnetic  field 
produced  in  the  E- region  could  have  an  amplitude  as  large  as  0.5  gamma, 
which  is  greater  than  most  natural  micropulsation  amplitudes  observed 
beneath  the  ionosphere . 

Preliminary  calculations  with  the  expressions  for  the  two  hm  wave 
modes  indicate  that  significant  hm  wave  amplitudes  can  be  produced  in  the 
ionosphere  above  the  larger  current  loops  considered  in  this  study.  This 
conclusion  is  independently  confirmed  by  Greifinger  (1972).  Thus,  we 
conclude  that  the  artificial  generation  of  micropulsations  is  feasible  and 
that  a  large  current  loop  should  be  capable  of  producing  observable  Pc  1 
micropulsations  at  very  considerable  distances  (greater  than  5000  km)  from 
the  source.  Much  further  study  is  required  to  establish  all  the  possibilities 
of  the  artificial  generation  process,  particularly  at  low  micropulsation 
frequencies  and  for  field- line  propagation.  However,  a  number  of  important 
uses  for  artificially-generated  micropulsations  are  immediately  apparent 
and  are  discussed  briefly. 
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I .  INTRODUCTION 


A.  DESCRIPTION  OF  MICROPULSATIONS 

Geomagnetic  micropulsations  are  small  approximately  periodic  fluctu¬ 
ations  of  the  earth's  magnetic  field  with  periods  in  the  range  0*2  -  600 

seconds  [see  reviews  by  Jacobs  (1970),  Troitskaye  (1967),  or  Campbell 

(1967)  ].  In  the  electromagnetic  spectrum  they  constitute  part  of  the  ULF 

(ultra-low-frequency)  band.  There  are  a  variety  of  naturally-occurring 

micropulsations  and  it  is  customary,  following  Jacobs  et  al  (1964) ,  to 

subdivide  the  micropulsations  into  two  main  classes  according  to  whether 

they  are  "of  a  regular  and  mainly  continuous  character"  (class  Pc)  or  are 

characterized  by  "an  irregular  pattern"  (class  Pi) .  These  two  main 

classes  of  micropulsations  are  then  further  subdivided  depending  upon 

their  period  as  shown  in  Table  1.1: 


Table  1.1 

Micropulsation  Nomenclature 


Notation 

Period  Range (sec,) 

Notation 

Period 

Range(sec . ) 

Pc  1 

0*2 

5 

Pi  1 

1 

40 

Pc  2 

5 

-  10 

Pi  2 

40 

-  150 

Pc  3 

10 

-  45 

Pc  4 

45 

-  150 

Pc  5 

150 

-  600 

The  division  into  subgroups  is  based  upon  the  physical  and  morphological 
properties  of  the  micropulsations  and  thus  classes  Pc  1  and  Pc  2,  for 
example,  represent  (or  are  intended  to  represent)  physically  different 
micropulsations.  This  system  of  classification  has  been  widely  adopted  by 
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micropulsation  workers.  It  should  be  noted,  however,  that  there  are 
certain  distinctive  varieties  of  micropulsation  activity  that  do  not 
fit  conveniently  into  the  standard  system  of  classification  and  have  their 
own  nomenclature.  Probably  the  most  important  of  these  special  varieties 
of  micropulsations  is  the  form  of  activity  known  as  'Intervals  of 
Pulsations  of  Diminishing  Period,'  or  IPDP's.  These  micropulsations  were 
first  defined  by  V.  A.  Troitskaya  (Troitskaya,  1961)  and  they  have  become 
the  subject  of  considerable  research  because  of  their  occurrence  during 
magnetospheric  substorms  [e.g.  Carpenter  et  al ,  1971,  and  references 
therein]. 

The  properties  of  the  various  classes  of  micropulsations  have  been 
investigated  rather  thoroughly  over  the  last  decade,  both  on  the  ground 
and  in  space.  It  is  now  clear  that  micropulsations  of  all  classes  have 
their  maximum  amplitudes  in  the  auroral  zones.  There  is  also  a  remarkable 
degree  of  correlation  between  micropulsations  recorded  at  geomagnet ically 
conjugate  points.  These  are  the  major  rimilarities .  Most  of  the  other 
properties  of  the  micropulsations  are  found  to  differ  markedly  from  one 
class  to  another. 

An  important  distinguishing  feature  of  the  Pc  1  class  of  micropulsations 
is  the  presence  of  a  fine  frequency  structure  that  is  easily  seen  on  their 
spectrograms  (i.e.  frequency-time  displays).  This  frequency  structure  is 
absent  from  all  other  classes  of  Pc  and  Pi  micropulsations.  Scattered 
structural  elements  are  sometimes  observed  during  IPDP  events,  but  they 
do  not  appear  to  be  a  primary  characteristic  of  the  micropulsations  as 
they  do  with  Pc  1,  Detailed  investigation  of  tho  frequency  structure  in 
Pc  1  micropulsations  resulted  in  the  present  theory  for  their  origin, 
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which  is  based  upon  a  particular  wave-particle  interaction  in  the  magneto¬ 
sphere  (e.g.  Cornwall,  1965).  it  is  not  IwJieveJ  that  wave-particle 
interactions  play  an  important  role  in  the  origin  of  any  of  the  other 
classes  of  Pc  or  Pi  micropulsations . 

Purely  electromagnetic  disturbances  at  micropulsation  frequencies 
do  not  propagate  in  the  ionospheric  or  magnetospheric  plasma.  The  waves 
that  do  propagate  involve  both  electromagnetic  fields  and  periodic  mass 
motions  and  they  are  strongly  influenced  by  the  ambient  geomagnetic  field  ; 
they  £  re  called  hydromagnetic  (hm)  or  magnetohydrodynamic  (mhd)  waves. 
There  is  an  extensive  literature  concerning  these  important  waves  and 
their  properties  will  not  be  discussed  here  in  great  detail.  The  render 
is  referred  to  the  micropulsation  review  text  by  Jacobs  (1970)  for  a 
relevant  discussion.  The  properties  of  micropulsations,  in  particular 
their  wide  occurrence  over  the  surface  of  the  earth  and  similarity  at 
conjugate  locations,  all  point  to  an  origin  as  hm  waves  in  the  ionosphere 
and  above.  Thus,  when  we  refer  to  the  artificial  generation  of  micro¬ 
pulsations  we  are,  strictly  speaking,  implying  the  generation  of  hm 
waves  with  micropulsation  frequencies  in  the  upper  atmosphere.  The 
surface  effects  of  these  hm  waves,  i.e.,  artificial  geomagnetic  micro¬ 
pulsations,  occur  when  the  hm  waves  reach  the  lower  part  of  the  ionospheric 
E-region  and  convert  to  electric  currents.  Magnetic  fluctuations  with  a 
superficial  similarity  to  natural  micropulsations  car.  be  produced  on  the 
surface  of  the  earth  by  various  means  (e.g.,  with  a  current  loop)  without 
intermediary  hm  waves.  These  disturbances  are  typically  extremely  localized 
and  are  physically  distinct  from  micropulsations.  We  will  not  use  the 
term  artificial  micropulsations  for  such  fluctuations. 
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B .  NEED  FOR  ARTIFICIAL  GENERATION 

Micropulsation  research  has  of  necessity  always  involved  the 
passive  acquisition  of  data.  Micropulsations  are  recorded  as  they  occur 
naturally  (and,  in  certain  cases,  infrequently)  and  the  records  are  then 
processed  to  obtain  information  about  the  micropulsations.  There  has 
never  been  a  wholly  controlled  method  of  generating  micropulsations. 
High-altitude  nuclear  explosions,  at  present  the  only  known  method  of 
artificial  generation  (e.g.  Green  et  al.,  1962;  Crook  et  al.,  1963; 

Bomke  et  al.,  1964;  Bowman  and  Mainstone,  1964;  Kovach  and  Ben-Menahem, 
1966),  produce  extremely  noisy  micropulsation  signals  in  addition  to 
having  other  obvious  disadvantages.  However,  insofar  as  micropulsation 
studies  are  concerned,  nuclear  explosions  have  two  very  desirable 
characteristics.  First,  the  source  of  the  micropulsations  is  localized 
and,  second,  the  generation  of  the  micropulsations  begins  at  a  precisely 
known  time.  These  characteristics  enable  propagation  times  to  be  calcu¬ 
lated  for  the  micropulsation  disturbances,  and  thus  Information  is 
obtained  about  the  hydromagnetic  waves  that  transmit  the  disturbances 
through  the  magnetosphere  and  ionosphere.  Micropulsation  propagation 
times  are  still  of  current  Interest  (see  Campbell  and  Thomberry,  1972), 
but  measurements  on  the  natural  signals  are  at  best  difficult  and  tedious. 

The  limited  usefulness  of  nuclear  explosions  for  the  generation  of 
micropulsations  suggests  that  a  fully  controlled  method  of  generation 
could  have  a  number  of  useful  applications.  (By  fully  controlled  we  mean 
that  the  frequency  of  the  micropulsation  source  as  well  as  its  location, 
timing  and  power  output  are  in  principle  all  at  the  disposal  of  the 
experimenter) . 


4 


This  report  is  concerned  with  the  possibility  of  using  a  very 
large  current  loop  on  the  ground  to  generate  hydromagnetic  waves  of 
micropulsation  frequencies  (i.e.  micropulsations)  in  the  overhead 
ionosphere.  There  are  many  advantages  associated  with  this  particular 
method  of  generation  compared  to  other  possible  methods.  Details  of 
the  current  loop  metnod  of  generation  will  be  discussed  in  the  next 
chapter.  At  this  stage  we  will  simply  list  some  of  the  more  important 
uses  for  an  artificial  source  of  micropulsations,  such  as  may  be  pro¬ 
vided  by  a  large  ground-based  current  loop: 

1)  Simple  and  direct  studies  of  hydromagnetic  wave  propagation 
in  the  magnetosphere  and  ionosphere  would  become  possible.  Precise 
studies  of  the  E-W  and  N-S  propagation  of  Pc  1  micropulsations  could  be 
conducted  for  the  first  time.  At  lower  frequencies  it  might  be  possible 
to  excite  resonances  either  along  an  entire  field  line  or  other  specific 
sections  of  a  field  line  (e.g.  between  the  lower  ionosphere  and  the 
Alfven  velocity  maximum  at  ~3000  km  altitude)  (e.g.  Jacobs  and  Watanabe, 
1362;  Fraser- Smith,  1967). 

2)  The  cyclotron  resonance  interaction  between  electrons  and  VLF 
waves,  and  between  protons  and  ULF  waves,  is  of  grea1:  current  interest 

in  magnetosphcrlc  studies.  Energetic  electrons  (protons)  may  be  precipi¬ 
tated  out  of  the  magnetosphere  by  VLF  (ULF)  wave  injection,  either  natural 
or  artificial,  as  a  result  of  this  interaction.  The  VLF  and  ULF  waves, 
in  turn,  may  be  amplified  through  their  interaction  with  the  energetic 
particles . 

At  present  there  is  little  direct  evidence  for  the  cyclotron  resonance 
interaction  and  its  importance  in  the  magnetosphere  is  largely  inferred  from 
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the  properties  of  certain  VLF  and  ULF  wave  phenomena  and  through  the 
observation  of  .limits  on  the  fluxes  of  trapped  charged  particles.  It 
is  of  interest  that  the  interaction  is  thought  to  be  essential  for  the 
generation  and  amplification  of  Pc  1  micropulsations  (e.g.  Cornwall,  1965; 
Liemohn,  1967) . 

Helliwell  et  al.  (1972)  recently  reported  some  of  the  most  direct 
evidence  yet  available  for  the  interaction:  the  observation  of  VLF 
whistler- induced  precipitation  of  energetic  electrons.  In  addition,  they 
noted  the  feasibility  of  conducting  controlled  experiments  with  a  high- 
powered  (100  kW)  VLF  transmitter  at  Siple  Station  in  the  Antarctic. 

A  large  ground-based  current  loop  would  have  a  unique  capability 
for  experiments  on  the  cyclotron  resonance  interaction  at  ULF  frequencies, 
i.e.  at  frequencies  much  less  than  in  the  planned  VLF  experiments.  In 
addition,  the  loop  could  precipitate  energetic  protons  out  of  the  magneto¬ 
sphere  and  into  the  lower  ionosphere  on  a  semi-controlled  basis .  These 
protons  would  be  useful  in  other  experiments.  Finally,  the  loop  could  be 
used  to  test  experimentally  the  reality  of  the  supposed  Pc  1  micropulsation 
generation  process . 

3)  Magnetotelluric  studies  of  the  crust  and  ipper  mantle  would  be 
facilitated.  These  studies  are  largely  dependent  upon  naturally-occurring 
low  frequency  electromagnetic  fluctuations  on  the  earth's  surface  and  there 
are  many  gaps  in  the  available  frequencies. 

4)  Provided  sufficient  power  is  available,  as  may  be  assumed  for  a 
ground-based  generator,  the  ionosphere  may  be  modified  either  by  direct 
heating  or  by  the  generation  of  currents  and  hydromagnetic  waves.  In 
other  words,  studies  of  the  ionosphere  may  be  possible.  At  frequencies 
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above  1  Hz  hydromagnetic  waves  are  strongly  absorbed  in  the  ionosphere 

Francis  and  Karplus,  1960;  the  location  of  the  maximum  of  absorption 
varies  with  frequency)  and  there  is  an  efficient  conversion  into  heat 
of  the  energy  deposited  in  the  ionosphere.  Unfortunately,  the  magnitude 
of  the  currents  induced  in  the  earth  by  the  current  loop  also  increases 
with  frequency.  These  currents  reduce  the  effectiveness  of  the  loop  and 
dissipate  energy;  they  will  be  discussed  in  greater  detail  in  Chapter  III. 

In  addition  to  the  above  uses  for  an  artificial,  controlled  source 
of  micropulsations,  there  are  less  tangible  benefits  to  be  expected  from 
a  first  attempt  to  conduct  an  experiment  upon  a  natural  system.  Much 
has  been  surmised  about  the  properties  of  the  ionosphere  and  magnetosphere 
at  micropulsation  frequencies  but  there  are  very  few  definite  experimental 
results.  Both  theoretical  and  experimental  work  would  be  stimulated  by  an 
attempt  to  generate  artificial  micropulsations,  and  there  may  be  new  dis¬ 
coveries,  For  example,  the  recent  successful  experiments  to  modify  the 
ionosphere  by  using  very  high  power  HF  radio  transmissions  (e.g.  Utlaut, 
1970;  Carlson  et  al.,  1972)  produced  some  results  that  were  completely 
unexpected  and  showed  that  the  ionosphere,  thought  by  many  to  be  well- 
understood  at  HF  frequencies,  still  had  some  novel  properties  that  could 
not  be  easily  explained. 

C.  PLAN  OF  THE  REPORT 

In  Chapter  II  an  outline  is  given  of  the  suggested  method  of 
generation  of  micropulsations  by  using  a  large  ground-based  current 
loop.  Chapters  III  and  IV  are  concerned  with  the  theory  for  two  of 
the  most  important  aspects  of  the  generation  process.  In  Chapter  III 
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the  effect  of  earth  currents  on  the  performance  of  the  loop  is  calculated 
and  in  Chapter  IV  expressions  are  derived  for  the  two  forms  of  Alfven 
wave  produced  by  the  loop  in  the  ionosphere.  Chapter  III  also  contains  a 
discussion  of  skin  depths  for  selected  materials  comprising  the  earth 's 
surface  and  crust.  Some  practical  details  of  the  loop  are  discussed  in 
Chapter  V  and  the  overall  conclusion  of  the  report  is  given  in  Chapter  VI. 
Three  appendices  contain  subsidiary  information  required  in  the  two 
theoretical  chapters. 
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II.  OUTLINE  OF  THE  CURRENT  LOOP  METHOD  FOR  GENERATING  MICROPULSATIONS 


A.  INEFFECTIVENESS  OF  CURRENT  LOOPS  IN  FREE  SPACE 

Loop  antennas  have  been  used  for  many  years  to  transmit  and 

receive  electromagnetic  signals .  In  most  cases  these  loops  operate  at 

frequencies  in  the  upper  kilohertz  or  megahertz  range  and  their  physical 

dimensions  are  small  relative  to  the  wavelengths  involved,  i.e,,  their 

radii  do  not  exceed  a  few  meters  .  Provided  the  loop  dimensions  are 

small  compared  to  the  wavelength  of  the  transmitted  signal,  a  loop 

antenna  used  as  a  transmitter  is  equivalent  to  a  magnetic  dipole  with  an 

oscillating  moment.  In  calculating  the  fields  radiated  from  such  a  dipole 

(e.g.,  Panofsky  and  Phillips,  1962)  the  wavelength  \  of  the  radiation  plays 

and  important  role  in  defining  the  radiation  and  induction  zones.  We 

will  use  the  term  "induction  zone"  to  mean  the  region  close  to  the 

dipole  (transmitter)  where  the  predominant  electric  and  magnetic  field 

-2  -  3 

components  vary  as  r  or  r  ,  with  r  indicating  the  distance  from  the 
dipole.  Another  term  commonly  used  for  the  induction  zone  is  "near¬ 
field  region".  The  "radiation  zone"  starts  at  some  distance  from  the 
dipole  and  extends  to  infinity.  In  this  zone  the  electric  and  magnetic 
field  components  vary  predominantly  as  r  1.  The  transition  region 
between  the  induction  and  radiation  zones  is  not  well-defined  but  it  is 
approximately  centered  on  the  radial  distance  r  =  \.  Thus  the  fields  are 
primarily  inductive  for  distances  less  than  a  wavelength  from  the  dipole 
and  for  distances  greater  than  a  wavelength  they  are  primarily  radiative. 
The  importance  of  the  wavelength  as  a  scale  factor  now  becomes  apparent, 
because  the  inductive  fields  do  not  contribute  to  the  radiated  energy 
from  the  dipole. 
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These  details  become  important  when  the  fields  generated  by  a 
large  ground-based  current  loop  are  considered.  For  a  1  Hz  frequency 

5 

the  wavelength  of  an  electromagnetic  wave  in  free  space  is  3  x  10  km, 
which  is  much  larger  than  the  radius  of  the  earth  (6370  km)  and  about 

4 

an  order  of  magnitude  larger  than  the  earth's  circumference  (4  x  10  km). 
Thus,  for  frequencies  in  the  micropulsation  range  the  direct  magnetic 
field  produced  by  a  current  loop  is  predominantly  inductive  at  any  point 
on  the  earth's  surface.  Close  to  the  loop,  i.e.,  within  a  few  hundred 
kilometers,  the  induction  magnetic  field  is  simply  the  standard  magnetic 
field  for  a  current  loop  with  an  e1^*  factor  to  allow  for  the  time 
variation  of  the  current  in  the  loop.  Expressions  for  the  magnetic 
field  produced  by  a  steady  current  in  a  conducting  loop  are  given  in 

_  3 

Appendix  A,  and  it  is  shown  that  the  field  drops  off  as  r  for  distances 
r  large  compared  to  the  loop  radius.  It  is  primarily  for  this  reason 
that  current  loops  in  free  space  have  been  considered  ineffective  as 
generators  of  magnetic  fluctuations  at  micropulsation  frequencies: 
there  is  no  significant  radiation  field  and  the  magnetic  field  that  is 
produced  by  the  loop  drops  off  too  rapidly  to  be  important  at  a  distance. 

B.  GENERATION  OF  HYDROMAGNETIC  WAVES  IN  THE  IONOSPHERE 

Although  the  magnetic  field  from  a  current  loop  drops  off  rapidly 
with  distance,  it  is  still  possible  for  a  large  current  loop  on  the 
ground  to  produce  a  magnetic  field  of  the  order  of  one  gamma  at  the 
height  of  the  ionospheric  E- region  (Details  of  such  loops  will  be  given 
in  Chapter  V) .  If  the  current  in  the  loop  is  varied  at  a  micropulsation 
frequency  the  E- region  will  experience  a  magnetic  disturbance  that  is 
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large  compared  to  most  natural  disturbances.  It  can  be  anticipated  that 
the  varying  magnetic  field  will  interact  with  the  E-region  plasma  to 
produce  hm  waves.  These  waves  then  propagate  away  from  the  interaction 
region,  possibly  to  large  distances,  either  within  the  ionosphere  or 
along  the  geomagnetic  field.  In  general,  artificial  geomagnetic  micro¬ 
pulsations  will  be  produced  on  the  ground  wherever  the  hm  waves  are 
strong  within  the  ionosphere,  either  through  ionospheric  ducting  (Pc  3 
frequencies;  see  Greifinger  &  Greifinger,  1968)  or  through  propagation 
along  the  geomagnetic  field  lines  (Pc  1  -  Pc  5  frequencies;  see  Jacobs, 
1970) . 

For  ionosp  eric  and  magnetospheric  plasmas,  and  micropulsation 
frequencies,  there  are  Just  two  modes  of  hm  wave  propagation  that  are 
of  importance.  These  two  modes  are  usually  called  "Alfvln"  and  "modified 
Alfven"  waves.  Both  modes  can  propagate  at  a  general  angle  to  the  ambient 
magnetic  field  but,  in  contrast  to  the  modified  Alfven  waves,  Alfven  waves 
tend  to  be  strongly  guided  by  the  field.  Also,  the  Alfven  waves  do  not 
propagate  iu  a  direction  perpendicular  to  the  ambient  magnetic  field,  as 
do  the  modified  Alfven  waves  (when  they  are  sometimes  referred  to  as 
magnetosonic  waves)  ,  The  velocity  of  Alfven  waves  along  the  magnetic 
field  is  the  well-known  Alfven  velocity  and  is  given  by 


where  Bq  is  the  steady  magnetic  field  and  p  the  plasma  density.  Two 
other  terms  commonly  used  for  the  Alfven  and  modified  Alfven  waves  are 
"slow  Alfven"  waves  and  "fast  Alfven"  waves,  respectively.  This  usage 
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reflects  a  difference  in  the  phase  velocities  of  the  two  waves.  However, 
provided  the  wave  frequencies  are  much  less  than  the  local  ion  gyro- 
frequency,  there  is  little  difference  between  the  velocities  of  the  fast 
and  slow  Alfven  waves  in  the  ionosphere  and  magnetosphere.  (The  use  of 
the  term  hydromagnetic  usually  implies  that  the  wave  frequencies  a  are 
much  less  than  the  ion  gyrofrequency  of  tl  -  plasma  propagating  the  waves) . 
We  expect  both  hm  modes  to  be  excited  by  the  current  loop  in  the  ionosphere, 
with  the  Alfven  mode  propagating  primarily  in  the  direction  of  the  geo¬ 
magnetic  field  and  the  modified  Alfven  mode*  propagating  primarily  within 
the  ionosphere  and  becoming  trapped  in  the  F-region  duct  (Greifinger  and 
Greifinger,  1968).  Expressions  for  the  two  modes,  as  excited  by  a 
ground-based  current  loop,  are  derived  in  Chapter  IV. 

It  should  be  emphasized  that  the  current  loop  does  not  itself 
radiate  artificial  micropulsations.  Instead,  energy  is  fed  into  the 
E-region  of  the  ionosphere  by  the  alternating  magnetic  field  of  the 
loop.  (The  magnetic  field  does  work  upon  the  ionosphere).  Hm  waves 
then  propagate  away  (radiate)  from  the  interaction  region  to  other 
parts  of  the  ionosphere  and  magnetosphere,  as  described  above.  There 
is  no  significant  radiation  of  micropulsation-frequency  electromagnetic 
waves  in  the  region  between  the  earth  and  the  ionosphere. 

It  will  be  quite  easy  experimentally  to  distinguish  the  artificial 
micropulsations  generated  by  the  current  loop  from  the  direct  fluctuating 
magnetic  field  of  the  loop.  The  direct  field  can  be  calculated  accurately 
and  at  moderate  distances  from  the  loop  the  magnetic  field  variation,  with 
artificial  micropulsations  present,  will  be  substantially  larger  than  the 
calculated  direct  field  variation.  At  large  distances  (2  1000  km;  the 
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exact  distance  will  depend  on  the  loop  magnetic  moment  and  the 
sensitivity  of  the  detection  equipment)  the  direct  field  will  be 
negligible  and  only  artificial  micropulsations  will  be  observed  together 
with  any  natural  micropulsations  that  might  be  present.  Conjugate  point 
observations  will  also  give  direct  verification  of  the  generation  process, 
because  the  direct  field  from  the  loop  will  probably  be  negligible  at 
the  conjugate  point.  Finally,  if  the  current  variations  in  the  loop 
are  coded,  the  artificial  micropulsations  (propagation  velocity  **>  V^) 
will  be  delayed  relative  to  the  direct  field  from  the  loop  (propagation 
velocity  ~c;  i  .e . ,  the  velocity  of  light  in  free  space).. 

The  critical  quantity  in  this  proposed  method  for  generating 
artificial  micropulsations  by  using  a  ground-based  current  loop  is  the 
magnitude  of  the  magnetic  field  that  can  be  produced  in  the  E- region. 

The  magnitude  of  the  magnetic  field  depends  on  only  two  factors:  the 
magnetic  moment  of  the  loop  and  the  conductivity  of  the  earth.  Of  these 
two  factors,  the  conductivity  of  the  earth  is  the  least  capable  of  being 
controlled.  The  effect  of  currents  induced  in  the  conducting  earth  is 
discussed  in  the  following  chapter.  There  appears  to  be  no  fundamental 
problem  in  constructing  current  loops  of  large  moment,  and  details  of  some 
suitable  loops  will  be  discussed  in  Chapter  V.  For  practical  reasons,  the 
requirement  for  a  large  magnetic  moment  restructs  consideration  to  current 
loops  that  are  coplanar  with  the  earth’s  surface. 
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III.  EFFECT  OF  INDUCED  EARTH  CURRENTS 


The  effect  of  induced  earth  currents  upon  the  micropulsation 
generating  capability  of  the  ground-based  current  loop  is  considered  in 
this  chapter.  The  induced  currents  are  important  because  their  magnetic 
field  opposes  the  inducing  magnetic  field  of  the  loop  and  the  total 
magnetic  field  is  everywhere  smaller  than  it  would  be  if  the  earth  was 
a  non-conductor.  Since  the  effectiveness  of  the  loop  as  a  micropulsation 
generator  depends  upon  the  total  field  it  produces  at  ionospheric  heights, 
it  can  be  seen  that  the  presence  of  the  conducting  earth  reduces  the 
capability  of  the  loop.  Fortunately,  the  earth's  surface  is  not  every¬ 
where  composed  of  materials  with  a  high  conductivity,  such  as  sea-water, 
and  some  important  material,  have  very  small  conductivities.  One  of  the 
major  objectives  of  this  chapter  is  to  find  how  strongly  the  direct  field 
of  the  loop  is  reduced  by  induced  earth  currents  for  a  wide  range  of 
conductivities  and  frequencies. 

The  chapter  starts  with  a  definition  of  skin  depth  and  a  brief 
description  of  its  variation  for  frequencies  in  the  range  10-3  -  103  Hz 
and  for  some  representative  earth-surface  conductivities.  A  complete 
analytic  solution  to  the  field  equations  is  then  presented  for  a  system 
consisting  of  a  current  loop  located  on  the  surface  of  a  single- layered 
earth.  The  plane  of  the  loop  is  taken  parallel  to  the  surface  of  the 
earth.  The  final  equation  for  the  total  field  is  solved  numerically 
and  graphs  of  the  variation  of  total  field  with  frequency  and  conductivity 
are  plotted. 

The  particular  theoretical  approach  used  in  this  chapter  is  basically 
the  same  as  that  introduced  by  Price  (1950)  and  later  extended  and  modified 

Preceding  page  blank 
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by  various  other  authors,  including  Price  (e.g.,  Gordon,  1951;  Bhattacharyya , 
1955;  Price,  1962),  There  is  also  a  very  considerable  literature  concerned 
with  related  problems  in  magnetotelluric  work  (e.g,,  Wait,  1962;  Mann,  1967; 
Ryu  et  al.,  1970).  It  should  be  pointed  out  that  there  is  a  very  different 
point  of  view  between  the  magnetotelluric  studies  and  the  work  reported 
here.  This  study  is  concerned  with  the  effects  produced  by  a  ground-based 
current  loop  in  the  ionosphere,  whereas  the  other  studies  are  motivated 
largely  by  the  possibility  of  sounding  the  structure  of  the  earth.  In  the 
latter  case  it  is  the  fields  on  the  surface  of  thr  enrth  (i.e.,  the 
measured  fields)  that  are  of  primary  interest  and  the  emphas i i  is  on  the 
behavior  of  the  surface  fields  when  the  earth  is  multi-layered  or  has 
other  structure,  such  as  a  discontinuity  in  the  horizontal  plane. 


A.  SKIN  DEPTH 

The  skin  depth  6  is  defined  by  the  equation 


6 


M0u,a 


(3.0) 


It  provides  a  convenient  measure  of  the  depth  of  penetration  of  an 
electromagnetic  field,  varying  with  angular  frequency  oj,  into  a  conductor 
of  conductivity  a  .  Equation  3,0  is  obtained  by  assuming  that  a  plane 
electromagnetic  wave  is  perpendicularly  incident  upon  the  surface  of  a 
plane  conductor  and  then  calculating  the  rate  at  which  the  electric  field 
falls  off  with  increasing  depth  in  the  conductor.  The  fall-off  is 
exponential  and  the  skin  depth  6  is  the  depth  at  hich  the  amplitude  of 
the  electric  field  has  fallen  to  l/e  of  its  value  just  inside  the 
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surface  of  the  conductor.  The  idealized  conditions  for  the  derivation 
of  Equation  3.0  do  not  always  apply  in  practical  situations  but  the  skin 
depth  defined  by  the  equation  still  provides  a  very  useful  indication  of 
the  depth  of  penetration  of  the  fields. 


Table  3.1 

Skin  Depths  for  Certain  Frequencies  and  Conductivities 


Frequency  (Hz) 

io"  4 

Conductivity  (mho/meter) 

IO-1  4 

-3 

10 

1590 

km 

50.3  km 

7.95 

km 

H1 

o 

1 

to 

503 

km 

15,9  km 

2.51 

km 

-1 

10 

159 

km 

5.03  km 

795 

m 

l 

50.3 

km 

1.59  km 

251 

m 

10 

15.9 

km 

503  m 

79.5 

m 

45 

7.50 

km 

237  m 

37.5 

m 

io2 

5,03 

km 

159  m 

25.1 

m 

103 

1.59 

km 

50.3  m 

7.95 

m 

Table  3.1  gives  the  skin  depths  for  a  number  of  frequencies  and 
conductivities  relevant  to  the  present  work.  Sea  water  is  easily  the 
best  conductor  (a  =  4  mho/m)  of  the  materials  occurring  in  large  quan¬ 
tities  over  the  surface  of  the  earth.  A  typical  surface  "earth"  conductivity 

-1  -2  , 

lies  somewhere  in  the  range  10  -  10  mho/m  with  the  conduction  caused 

largely  by  the  presence  of  saline  water.  For  certain  dry  rock  and  sand 
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—4 

formations  the  conductivity  may  fall  below  10  mho/m.  Pure  ice  and 
pure  water  have  very  low  conductivities  (of  the  order  of  10  -  10  mho/m) 

but  the  conductivities  observed  in  nature  for  ice  and  "pure"  water  vary 
widely  depending  upon  the  content  of  impurities,  especially  salt. 

Table  3.2,  derived  from  data  published  by  Evans  (1965),  lists  some  con¬ 
ductivities  for  ice,  snow  and  water. 


Table  o»2 

Conductivities  of  Ice,  Snow  and  Water 


Material 

Conductivity  (mho/meter) 

Pure  Ice  (-10®C) 

IQ"7 

Glacial  Ice,  Alberta  (0®C) 

5  x  10  7 

Sea  lex  (0.5%  salinity;  -7°C) 

CO 

1 

o 

rH 

Soft  New  Snow  (-10°C) 

10-9 

Compact  Wet  Snow  (0°C) 

10-6 

Pure  Water 

4  x  10"6 

Distilled  Water 

2  x  10"4 

Public  Water  Supply 

10'3  -  5  x  10“2 

Sea  Water 

1-5 

The  earth's  crust  is  generally  considered  to  be  shout  three  miles 
thick  (~  5  km)  under  the  sediments  of  the  ocean  floors  and  perhaps  twenty 
miles  thick  (~30krd),  on  the  average,  under  the  continents.  A  layer  of 
sedimentary  rocks  with  a  thickness  varying  from  0  -  5  km  forms  the  upper 
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part  of  the  continental  crust.  The  sedimentary  rocks  are  comparatively 

highly  conducting,  with  conductivities  in  the  range  10  -  10  mho/m 

as  noted  earlier.  The  magnitude  and  variation  of  the  conductivity  below 
the  sedimentary  layer  has  been,  and  still  is,  the  subject  of  considerable 
debate.  A  good  introduction  is  given  by  Mott  and  Biggs  (1963).  We  will 
only  comment  that  the  comparatively  cool  granitic  (upper)  and  basaltic 
(lower)  materials  comprising  the  major  part  of  the  crust  have  very  low 
conductivities,  possibly  as  low  as  10  7  mho/m.  The  Mohorovicic  discon¬ 
tinuity  separates  the  crust  from  the  aathenosphere,  or  upper  part  of  the 
mantle.  Below  this  discontinuity  there  is  an  increase  in  conductivity 
and  it  may  reach  values  in  excess  of  10  ^  mho/m. 

A  conductivity  of  10-4  ;aho/m  appears  a  reasonable  estimate  for  the 
crust  of  the  earth  at  continental  locations  where  the  sedimentary  layers 
are  either  thin  or  of  low  conductivity.  Table  3.1  shows  that  the  varying 
magnetic  field  of  the  loop  will  easily  penetrate  the  crust  if  the  frequen¬ 
cies  are  about  1  Hz  or  less.  For  ~1  Hz  fluctuations  (1  ,e.,  Pc  1  micro¬ 
pulsations)  it  is  probably  still  a  good  first  order  approximation  that 
the  earth  is  a  single  layer  of  conductivity  10  4  mho/m,  but  for  frequencies 
less  than  about  0.2  Hz  a  two  layer  model  would  be  more  appropriate.  The 
two  layer  profile  of  Cantwell  and  Madden  (1960)  for  Massachusetts,  con¬ 
sisting  of  an  upper  layer  70  km  thick  and  of  conductivity  1.2  x  10  4  mho/m 
above  a  basement  layer  of  conductivity  less  than  1.2  x  10  *  mho/m,  may  then 
provide  a  good  starting  point. 
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B.  FIELD  INDUCED  BY  A  CURRENT  LOOP 

The  field  equations  are  now  solved  for  a  current  loop  source 
located  above  a  single- layered  conducting  earth.  The  earth  is  represented 
by  a  semi- infinite  homogeneous  medium  of  uniform  conductivity  o  occupying 
the  half-space  z  <  0  of  a  cylindrical  coordinate  system  (r,  <t> ,  z)  .  The 
semi-infinite  region  above  the  earth  (z  >  0)  is  assumed  to  be  free  space 
(conductivity  a  =  0,  permittivity  CQ,  permeability  |iQ).  It  is  also 
assumed  that  the  earth  is  non-magnetic ,  i.e.,  that  the  permeability  is 
everywhere  the  same  as  that  of  free  space.  Initially,  generality  is 
maintained  by  specifying  only  that  there  is  some  field  source  at  a 
height  h  above  the  surface  of  the  earth  (z  =  o) ,  that  the  z-axis  passes 
through  the  source,  and  that  there  is  rotational  symmetry  of  the  electro¬ 
magnetic  fields  about  the  z-axis.  In  most  previous  work  the  field  source 
is  assumed  to  be  a  system  of  ionospheric  currents  and  h  is  the  height  of 
the  lower  ionosphere.  In  the  present  work  h  will  ultimately  represent 
the  height  of  the  current  loop  above  the  ground,  and  the  situation  of 
most  interest  occurs  when  h  =  o  and  the  loop  is  actually  located  on  the 
ground . 


1 .  Basic  Equations 

The  following  are  the  expressions  for  grad,  div,  curl  and  divgrad 
in  cylindrical  coordinates: 


grad 


Y  = 


3Y  ~  1 

3F  r 


3Y 


A 

Z 


(3,1) 
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(3.2) 


div  A 


i  a<rV 

"  r  5r 


curl  A  = 


•  3A  dA  i 
r  z 

w '  'srJ 


L- 


ra<rV 


dr 


dA  n 
r  1 

&0 


z 


(3.3) 


divgrad  Y  = 


•> 

3? 


(r 


dY 


a2Y  52Y 

302  +  dz2 


(3.4) 


where  r  #  2  are  unit  vectors  in  the  appropriate  coordinate  directions . 

/*w^  <**•«/  ^ 

It  is  permissible  at  micropulsation  frequencies  to  ignore  the  dis¬ 
placement  current  density  D  and  no  generality  is  lost  by  assuming  that  the 
initial  space-charge  distribution  is  everywhere  zero  (Price,  1950).  The 
field  equations,  in  MKS  units,  are  then 

div  D  =  0 

div  B  =  0 

curl  H  =  J 
curl  E  =  -6 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


For  the  particular  source  considered  in  this  problem  there  is  rotational 
symmetry  of  the  electromagnetic  fields  and  it  is  easy  to  show  that  the 
components  B0,  Er  and  Ez  must  all  be  zero.  The  magnetic  and  electric  fields 
thus  consist  of  the  three  components  as  follows 
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B  =  <Br,  0,  Bt) 
E  =  <0,  E#,  0) 


The  coordinates  and  orientation  of  the  field  components  are  shown  in 
Figure  3.1. 

Solutions  are  now  derived  for  the  non-zero  field  components 

B  ,  B  and  E , .  (Because  there  is  only  one  component  of  electric  field 
r  z  0 

the  subscript  0  will  be  dropped  from  E^) .  The  time  variation  of  the 
electric  and  magnetic  fields  is  taken  to  be  sinusoidal  of  angular 
frequency  id  and 


where  the  amplitudes  E  and  B  are  functions  of  position  only. 

Initially,  consider  the  fields  within  the  conducting  medium.  Equations 
(3.7)  and  (3.8)  give 


E 


(3.9) 


n  1  aE 

r  Jto  "5z 


(3.10) 


B 

z 


_i  /  3e 
Jw  ("Sr 


(3.11) 
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Combining  these  equations  gives 


32e  32e  1  ar  n 


e  a  e  1  ar  fi  .  i  „ 

?  +  ^7  +  ■'*  "  L7  *VJ  E  =  0 


(3.12) 


where  ujh^ct  is  twice  the  inverse  square  of  the  skin  depth. 

Solutions  to  (3.12)  are  obtained  by  the  method  of  separation  of 


variables .  Let 


E  =  R(r)  Z(z) 


(3.13) 


Substituting  (3.13)  in  (3.12)  gives 


1  d^Z 
1  dz2 


4  .  / 1  d  R  1  dR  l\ 

**>  )  (*  ^7  57  ®  ‘  7) = 


=  0  . 


Because  the  quantities  in  parentheses  are  independent  it  follows  that 


A  ,1  *  ,  r,2  n  . 

7*55*'  ‘  ~sl  R  * 

dr  *-  r  J 


(3.14) 


^-|  -  (X2  +  JuHqC)  Z 

dz 


(3.15) 


where  -X  is  the  separation  constant  and  where  the  gj^ct  is  twice  the 
inverse  square  of  the  skin  depth. 

Equation  (3.14)  may  be  written  in  the  form 


d(Xr)'1 


(Xr)  d(Xr) 


*  (‘  -  ^7)  "  =  ° 


(3.16) 
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which  is  a  form  of  Bessel 's  equation  with  solutions 


R(r)  =  Jx(Xr)  (3.17) 

i.e.,  Bessel’s  functions  of  the  first  kind,  of  order  1. 

There  are  two  solutions  to  equation  (3.15).  For  0<  z<h,  i.e.,  in 
free  space,  ct  =  o  and 


Z  =  C/Z  +  C2e"Xz  .  (3 

For  z  <  0,  i.e.,  within  the  earth,  it  is  only  possible  to  specify  solutions 
Z(z,\)  that  vanish  at  infinity. 

Thus, 


=  Z(z,X)  J^(Xr)  z  <  o 

E+  =  +  C2e  ^z)  J^Xr)  0  <  z  <h  .  (3.19) 

The  standard  boundary  conditions  for  the  interface  at  z  =  o  reduce 

to  the  conditions  that  E,  B  and  B  should  all  be  continuous  across  the 

r  z 

boundary.  Application  of  these  conditions  gives 

"  Ci  +  C2  (3.20) 

■  "  C2}  *  (3.21) 


Z(0,  X) 


Ptel 
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By  using  (3.20)  and  (o.21)  and  the  relationship 


J1  (x)  ®  J0(x)  "  ^  Ji(x) 


two  equations  for  the  magnetic  field  components  in  the  region  h>z>o 
are  obtained,  i.e., 


Br  =  +  <Cie>'2  -  J1  <X-'> 


(3.22) 


i;  <cie>'2  ♦  v"*“>  vXr) 


(3.23) 


Also,  since 


J0(x) 


Equations  (3.22)  and  (3.23)  can  be  written  in  the  alternate  form 


,  3  r  Jo(  ^  xz  „  -  \z.  "i 

Jr  *  ITr  [~  ~  <V  -  V  >J 


(3.24) 


r  J0Ur)  \Z  „  -Xzl 

[-  ~ —  <Cle  C2e  >J  ’ 


B  =  -5- 

z  oz 


(3.25) 


In  free  space  curl  £  =  0  and  a  magnetic  scalar  potential  Q  may  be 
defined  by  the  equation 


B  =  -  grad  n 

U 


(3.26) 
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Equations  (3.24)  and  (3.25)  show  that  in  this  problem  the  magnetic 


scalar  potential  in  the  region  0<z  <h  is  of  the  form 


a 


J  (Xr) 
0 

JojMq 


V'*2) 


(3.27) 


AZ 

The  part  of  Cl  that  involves  e  increases  with  z  (in  the  indicated  region) 

and  may  be  identified  with  the  inducing  field,  while  the  part  involving 
-  Xz 

e  represents  the  field  due  to  induced  currents  in  the  conducting 
region  (e.g.,  Price,  1950,  1962;  Battacharyya,  1955;  Weaver,  1964).  The 
induced  field  is  finally  found  by  expressing  C2  and  Z  (z,  X)  in  terms  of 
C^,  which  is  assumed  to  be  known.  The  present  problem  is  not  directly 
concerned  with  the  induced  field  within  the  conductor.  However,  it  is 
still  necessary  to  consider  the  Z(z,  X)  solutions  in  3ome  detail  in 
order  to  relate  C2  to  .  The  appropriate  form  of  for  a  current  loop 
source  must  also  be  found  or,  equivalently,  the  inducing  potential  must 
be  expressed  in  the  form  of  a  sum  (or  integral)  of  terms  of  the  type 
Ce^Z  JQ(Xr) . 


2 .  Effect  of  the  Earth 

The  solutions  Z(z,X)  to  Equation  (3.15)  will  now  be  considered  in 
greater  detail.  The  equation  may  be  written 


4  -  e2z  =  o 

dz2 


(3.28) 
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where 
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Thus,  outside  the  earth  (0<z<h)  the  magnetic  scalar  potential  of  the 
induced  magnetic  field  is  given  by 


a 


i 


(9-\) 

75+X) 


C1  e  J0(kr) 


(3.32) 


3 .  Current  Loop  as  Source 

The  current  loop  dimensions  considered  in  this  report  are  mostly 
small  compared  to  the  earth-ionosphere  distance  and  a  dipole  approach 
is  adequate  for  calculation  of  the  electromagnetic  fields.  A  plane 
current  loop  carrying  current  I  and  of  infinitesimal  are-  dS  has  a 
magnetic  dipole  moment  M,  where 


In  free  space  the  magnetic  scalar  potential  for  such  a  dipole  is  given  by 
n  _  M  (z  -  h) 

d  4ir  r  TTf  .■.2.3/2  (3.34) 

[r  +(z-h)  ]  ' 

where  the  dipole  is  located  at  z  =  +h  on  the  z  axis  of  a  cylindrical 
system  of  coordinates  and  has  its  axis  aligned  along  the  (positive)  z 
direction.  Equation  (3.34)  can  be  written 


where  p  =  [r2  +  (z  -  h)2]1/2  is  the  radial  distance  from  the  dipole  to 
the  field  point.  The  following  formula  (Sneddon,  1955) 

00 

f  JQ(Xr)  e^<Z"h)dX  =  i  (z<h)  (3.36) 

0 

may  be  used  to  write  (3.35)  in  the  form 


JQ(Xr)  e 


X(z-h) 


\d\ 


(z<  h) 


For  a  dipole  at  the  origin,  and  z>0, 


B 

r 


B 

z 


H-M  “  1  o 

~  f  J.(Xr)e"Xz  X2dX 
4jt  J  1 
0 


(3.37) 


(3.38) 


The  negative  sign  in  equation  (3.37)  comes  from  the  assumption 
that  the  magnetic  dipole  moment  is  directed  in  the  positive  z  direction. 
The  dipole  moment  is  actually  oscillatory  (I  =  I0e^wt)  and  the  sign  is 
not  really  significant. 

The  potential  given  by  equation  (3.37)  satisfies  the  requirements 
previously  discussed  for  a  suitable  inducing  potential.  Combining  (3.37) 
and  (3.32)  the  following  equation  for  the  induced  field  in  the  region 
z  >0  is  obtained 
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(3.39) 


M  00 

Of  8-X  -\(z+h)  ... 

h  -  TR  J  ' 571  e  Jo(Xr)  XdX 

o 

where  MQ  =  IQdS .  Thus,  for  the  dipole  (loop)  at  the  origin  of 
coordinates,  the  induced  field  in  the  region  z  >0  is  given  by 

n  Mo  r\x  -Xz  ■  ,  , 

ni  =  Tn  J  e  J0(Xr>  ^  •  (3.40) 

0 

Equations  (3,39)  and  (3.40)  represent  the  formal  solution  to  our  problem. 


4 .  Numerical  Solution  of  the  Induced  Field  Equation 

We  now  concentrate  upon  Equation  (3.40)  for  the  induced  field 
produced  by  an  infinitesimal  current  loop  placed  upon  the  earth's  surface. 
The  integral  in  this  equation  was  solved  analytically  by  Gordon  (1951), 
but  unfortunately  the  solutions  so  obtained  are  unwieldy  and  limited  in 
their  application.  We  chose  to  evaluate  the  integral  numerically  by  using 
one  of  the  integration  routines  available  for  Stanford's  IBM  360/67 
computer.  First,  however,  it  was  necessary  to  derive  a  suitable  expression 
for  the  real  part  of  the  induced  field  equations.  Equations  (3.40)  and 
(3.26)  give  the  following  two  equations  for  the  induced  magnetic  field 
components: 


B 

r 


B 

z 


(3.41) 


Note  that  for  r  =  0  (i  .e .  on  the  z  axis)  we  have  ^(0)  =  0  and  JQ(0)  =  1. 
Thus,  Bf  =  0  as  required  and 

||  00 

....  0  f  0-X  -Xz  , 

n  (axis)  =  —  J  ^  e  XdX 

0 

(3.42) 

00 

«0  C  9-X  -Xz  .2  . 

\  <9Xls)  =  5/  wx  9  x  dX  • 

The  solution  of  the  equation  for  B  (axis)  will  enable  us  to  compare 

z 

the  direct  field  from  the  loop  with  the  induced  field  in  the  most  important 
region  of  the  ionosphere  for  the  production  of  hydromagnetic  waves,  i.e., 
in  the  E- region  directly  above  the  loop. 

Let 

0  =  a  +  jb  =  (X2  +  Ju^a)1/2  (3.43) 


where  both  a  and  b  are  real.  Solving  for  a  and  b,  we  obtain 

-1/2  2  4221/2  1/2 

a  =  2  [X2  4  (X4  +  )1/2] 

-^/2  2  4  2  2  2  l/2  ^2 

b  =  2  [-X  +  (X4  +  un0  aT'  j 


(3.44) 


The  real  part  of  equations  (3.41)  are  separated  out  by  writing 


8-X  _  (a-X)  +  jb 

0+X  ia+A.)  +  jb 

a2vb2-X2  2Xb 

(a*X)2+b2  +  (a+X)2+b2 
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Thus,  we  obtain  the  final  expressions  for  B  and  B  : 

r  z 


B  (real) 
r 


B  (real) 

z 


40M0  f[  a2+b2-\2  1 
4*  i  (a+X)2+b2 


e  Jg(Xr)X2  dX 


(3.45) 


^0M0  (Ta2+b2-X2  I 

4*  {  L(a+X)2+b2J 


2-72  dX 


with  the  auxiliary  quantities  a  and  b  defined  by  equations  (3.44)  . 

Note  that  as  a  -»  0  both  B  (real)  and  B  (real)  also  tend  to  zero, 

and  as  G  -t  00  we  obtain  simple  dipole  expressions.  J.u  other  words,  the 

induced  fields  given  by  equations  (3.45)  vary  in  the  required  manner. 

The  integral  in  the  equation  for  B  (real)  with  r  =  0  (i.e,,  axial 

z 

p 

values)  was  integrated  numerically  for  frequencies  in  the  range  10  ~  100  Hz 

-5 

and  conductivities  in  the  range  10  -  4  mho/m.  Other  data  used  in  the 

11  2 

computation  were  z  =  100  km  and  MQ  =  jt.10  amp.  m  .  These  latter  values 
correspond  to  a  field  point  at  E-region  height  and  a  loop  of  10  km  radius 
carrying  a  current  of  1000  amp.  With  no  contribution  from  the  earth  the 
magnitude  of  the  Bz  field  at  the  field  point  would  be  62.83  my  (see 
Appendix  A) . 

Figure  3.2  summarizes  the  results  obtained  from  the  numerical  inte¬ 
gration.  The  total  B^  field  is  plotted  in  the  figure  and  each  curve  shows 
the  variation  of  the  total  field  with  frequency  for  a  particular  value 
of  the  earth  conductivity.  The  dashed  line  indicates  the  "ideal"  or  zero 
earth  conductivity  magnitude  of  the  field.  It  cen  be  seen  that  for  a 
good  conductor  such  as  sea-water  (a  =  4  mho/m)  the  total  field  is  very 
much  less  than  the  ideal  value;  at  1  Hz  the  total  field  is  less  than 
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one  percent  of  the  zero  conductivity  field.  However,  for  conductivities 

typical  of  the  solid  earth  the  reduction  of  the  dipole  field  by  the 

-4 

induced  field  is  not  nearly  so  marked.  For  a  conductivity  of  10  mho/m, 
which  we  consider  typical  of  some  regions  of  the  earth,  the  total  field 
is  51.4  my  for  a  frequency  of  1  Hz.  This  value  is  only  about  18  percent 
smaller  than  the  zero  conductivity  field. 

At  higher  frequencies,  and  especially  for  the  greater  conductivities, 
the  curves  in  Figure  3.2  become  straight  lines.  In  fact,  as  shown  in 
Appendix  B,  when  5s  ^  skin  depth  the  total  field  becomes  proportional  to 

-1/2 

(uicr)  .  When  the  conditions  for  this  approximation  apply,  and  for 

the  data  used  to  prepare  Figure  3.2,  the  B field  may  be  written 

2  37 

B  (total;  =  - mY  (z  >  skin  depth)  .  (3.46' 

(coa)1/2 

This  equation,  or  the  more  explicit  form  in  Appendix  B,  may  be 

useful  under  some  circumstances.  However,  it  is  only  very  approximate  for 

-4 

the  frequencies  (~  1  Hz)  and  conductivities  (~  10  mho/m)  that  are  of 
particular  interest  in  this  report. 

C .  CONCLUSION 

The  results  presented  in  this  chapter  enable  us  to  reach  an  important 
conclusion  about  the  probable  effect  of  induced  earth  currents  on  the 
micj  ^pulsation  generating  capability  of  the  ground-based  current  loop. 
Provided  some  care  is  taken  in  choosing  the  location  of  the  loop,  it 
appears  likely  that  the  magnetic  field  produced  by  the  earth  currents  will 
remain  considerably  smaller  than  the  direct  field  of  the  loop.  For  a 


34 


TOTAL  MAGNETIC  FIELD  Bj  (mllllgamma) 


lo',o- 


FREQUENCY  (Hz) 


Figure  3,2 


Frequency  variation  of  the  maximum  axial  magnetic 
field  produced  by  the  10  km  loop  (maximum  current 
1000  amps)  at  an  altitude  of  100  km.  The  parameter 
a  (mho/m)  represents  the  conductivity  of  the  earth, 
considered  to  be  a  single  layer,  and  the  dashed  line 
indicates  the  field  for  a  non-conducting  earth. 


frequency  of  I  Hz,  if  the  loop  is  located  on  a  reasonably  non-conducting 
surface  material  (a  «  10~3  -  10~4  mho/m  in  the  first  1-5  km  of  crust),  it 
is  probable  that  the  magnetic  field  of  the  loop  at  E-region  heights  will 
not  be  reduced  by  more  than  20  percent  from  its  zero  earth  conductivity 
value . 

A  more  extended  calculation  of  the  effect  of  induced  earth  currents 

is  desirable  but  it  is  unlikely  to  affect  the  above  conclusion  to  any 

significant  extent.  A  more  realistic  model  of  the  earth’s  structure  would 

primarily  require  the  addition  of  a  deep  layer  of  higher  conductivity  than 
—4 

10  mho/m.  The  additioi  of  this  deep  layer  would  produce  ics  greatest 

effect  at  the  low  frequencies  (less  than  1  Hz)  where  the  skin  depth  is 

comparatively  large.  However,  the  curves  in  Figure  3.2  show  that  the 

magnetic  field  of  the  loop  at  low  frequencies  is  negligibly  affected  by 

-4 

the  upper  layer  of  conductivity  10  mho/m  and  the  depth  of  the  more- 
conducting  layer  will  necessarily  minimize  the  field  of  its  induced 
currents  at  ionospheric  heights  (the  field  from  the  induced  currents 
probably  drops  off  approximately  as  the  inverse  cube  of  the  distance  to 
the  field  point) . 
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IV.  GENERATION  OF  ALFVEN  WAVES  IN  THE  IONOSPHERE 


In  this  chapter  we  derive  expressions  for  the  Alfven  wavss  excited 
in  the  ionosphere  by  a  current  loop  on  the  ground.  In  order  to  use 
Fourier  techniques  we  utilize  the  model  shown  in  Figure  4.1  consisting  of 
two  adjoining  semi-infinite  homogeneous  regions  permeated  by  a  constant 
magnetic  field,  BQ.  The  lower  region  is  free  space  and  contains  the 
current  loop,  while  the  upper  region  is  a  plasma  and  supports  the  propa¬ 
gation  of  Alfven  waves. 


A .  COORDINATE  SYSTEMS 

Three  separate  coordinate  systems,  shown  in  Figure  4.2,  are  used 

to  study  the  waves  associated  with  the  problem.  The  (x,y,z)  system  is 

chosen  so  that  the  z-axis  points  in  the  vertical  direction  and  B  lies 

in  the  (x  ,  z  )  plane.  The  (x',y',z')  system  is  chosen  so  that  z  and  z’ 

coincide  and  so  that  Jk  lies  in  the  (x',z')  plane.  A  third  coordinate 

system,  (x" ,  y" ,  z")  is  chosen  so  that  B.  points  in  the  z"  direction  and 

'HJ 

k  lies  in  the  (x",z")  plane.  B  and  k  can  be  written  in  these  systems 

^  rsj 

as 


B„  =  B.(sin  cp  i  4  cos  cpi  ) 
0  Y  ~x  Y~z 


=  v1^ +  ra^ +  ni;) 


(4.1) 


B-i" 

0~z 


k  =  k  i  +  k  i  +  k  i  =  k(sin  9  i'  4-  cos  9  i') 
-  x~x  y~y  z~z  ~x  ~z 

=  k( sin  Y  i"  +  cos  Y  i") 

~x  z 


(4.2) 


where  i  ,  i  .  etc.  are  unit  vectors. 
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The  primed  and  unprimed  coordinate  systems  are  related  by  the 
transformation 


r  = 


k 

_y 

k. 


2  2 

where  kA  =  (kx  +  ky  ) 


1/2 


(4.3) 


Furthermore,  it  can  be  shown  that  the  primed  and  double-primed  coordinate 
systems  are  related  by  the  transformation 


r '  = 


sin  Q(l-t  )  -  to  cos  9  m  cos  0 
sin  Y  sin  * 

n  sin  9-1  cos  6 


m 


ctn  Y 


m 


cos  0(l-n  )-(zicos9 
sin  Y 


sin  Y 

m  sin  9 
sin  Y 


...  (4.4) 


By  reference  to  Figure  4,2,  it  can  also  be  demonstrated  that 


l  =  —  sin  cp 


m  =  r—1  sin  cp 


n  =  cos  cp 


sin  0  =  kA/k 

cos  Y  »  l  sin  0  +  n  cos  0  = 


(4.5) 


(4.6) 


k  sin  cp  sin  0 

— -  +  cos  0  cos  cp  (4.7) 
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B .  FREE  SPACE  WAVES  DUE  TO  CURRENT  LOOP 

The  governing  equations  in  free  space  are 


7  x  l  =  -JwHqH 
7  x  £  =  jw«o£  +  £ 

Their  Fourier  transform  is  given  by 

fe  *  4  -  “>*<& 

-J£  *  4  -  jo,eo4  ♦  4 

or,  upon  elimination  of  E  by 

~k  ’ 


O^k 


where  k  =  — 
0  2 
c 


(4.8) 

(4.9) 


(4.10) 

(4.11) 

(4.12) 

(4.13) 


Taking  the  scalar  product  of  this  equation  with  k  gives 

*4  =  j«m0  &-4A02 

and,  finally, 


E.  =  - 
~k 


JaVko  4  - 

ko2<k2  -  ko2> 


In  Appendix  C  it  is  shown  that 

4:  =  ~3  x  &  exP  ("Jk  a) 
8n  z 


(4.14) 


(4.15) 


(4.16) 


where  &  is  the  vector  magnetic  moment  of  the  current  loop.  Therefore  E 

l 

in  equation  (4.15)  becomes 


<4*0  JS,  x  & 

4  -  -  “3  -T— 2  e*P  [-Jkzal 
8*  k  -ko 


(4.17) 
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\ 


Taking  the  inverse  transform  with  respect  to  z  gives 


&kx.ky.z) 


“  X  5J, 

3  2  ~-2  exP  [-Jkz(z+a)  ]dkz 

8  it  k  -kQ 


(4.18) 


Performing  the  usual  contour  integration  gives 


E(k  k  z)  = 
~  x  y 


/  2  '  2 

where  k^  =  /kQ  -kA 


Jwp0  &  x  h  exP  i-Jkz0(z+a)  3 


„  2 
8jt  k 


(4.19) 


(4.20) 


For  future  reference,  note  that  this  may  be  written  in  the  primed 


coordinate  system  as 


E  =  E'  exp  f-jk  _(z+a)]  i' 
~  yi  F  1  J  zO  J 


(4.21) 


where 


-JkJ.wti0M 


(4.22) 
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c .  ALFVEN  WAVES 

The  theory  for  Alfven  waves  is  well  documented  and  only  the  results 
will  be  quoted  here  (Clemmow  and  Dougherty,  1969).  The  fast  wave  obeys 
the  dispersion  relation 

(4.23) 


where 


L1  " 


1  + 


fe) 


B„ 


VA  “ 


(4.24) 


VA  is  the  Alfven  velocity,  mg  and  m.  are  the  electron  and  ion  mass,  and 

N  is  the  equilibrium  number  density.  Referring  to  Figure  4.2,  the  electric 

field  is  linearly  polarized  in  the  y"  direction  i  e  E  -  E"  i" 

’  ’ '  ~  ~  y  ~y 

The  slow  wave  obeys  the  dispersion  relation 


where 


(4.25) 


(4.26) 


and  w  is  the  plasma  frequency. 

The  H-field  is  linearly  polarized  in  the  y"-direction  and  the  E-field 
is  elliptically  polarized  in  the  (x,z)  plane  such  that 


E" 

X 


k  cos  ¥ 


weoLi 


H" 

y 


E" 

z 


-k  sin  Y 


weoLo 


(4.27) 
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D.  SOLUTION  OF  THE  BOUNDARY  VALUE  PROBLEM 

When  each  plane  wave  generated  by  the  current  loop  strikes  the 
boundary  with  the  plasma,  two  plane  waves  are  reflected  b^ck  into  free 
space  and  fast  and  slow  Alfven  waves  are  transmitted  into  the  plasma. 
In  this  section  the  amplitudes  of  the  Alfven  waves  are  calculated. 

1 .  Fast  Alfven  Wave 

We  now  consider  the  fast  and  slow  Alfven  waves  and  designate 
quantities  relative  to  them  with  subscripts  a  and  b,  respectively.  The 
fast  wave  has  only  one  component  Ey^'  in  the  double  primed  system,  and 

therefore,  using  Eq .  4.4,  its  electric  field  components  in  the  primed 
system  are  given  by 


E' 

xa 


m  c°s  8  E  „ 

sin  y  ya 

a 

n  sin  0  -  l,  cos  0 

_  a  a 


sin  Y 


E' 

za 


m  sin  0 
_ a 

sin  Y 

a 


(4.28) 

(4.29) 

(4.30) 
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2. 


Slow  Alfv£n  Wave 


Now  consider  the  slow  wave.  Since  there  Is  only  one  component 
of  H^in  the  double-primed  coordinate  system,  H^,  we  find  from  Eq .  (4.4) 
that  the  H-field  in  the  primed  system  is  given  by 


m  cos  0 

H«  =  - 2  H  " 

xb  yb 

sin  Yw  3 
b 


(4.35) 


n  sin  9.-1  cos  6,. 

H.  - - b - b  H  " 

yb  sin  »  yb 

D 


(4  .36) 


m  sin  6 

H*  - - -  H  " 

zb  .  w  yb 

sin  Y.  3 

D 


(4.37) 


Likewise,  upon  substituting  Eqs .  (4.27)  into  Eq.  (4.4),  we  obtain 


E'  = 
xb 


k  cos  Y.  sin 
_ b _ b 

weoLi  Sln  Yb 


k  sin  Y» 


(4.38) 


W  JJ  M 

k  sin  Y.  ^ 
b 


(4.39) 


k  cos  Y.  cos  nu)u 

_ b _ b 

weoLo  sin  Yb  k  sin  Yb 


(4.40) 
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3 .  The  Incident  Wave 

The  incident  wave  has  only  one  component,  E^,  *hich  lies  in 
the  y'  direction.  From  Eq .  (4.8)  its  H-vector  has  the  components 


- - cos  0  E 

xl  oj|an  I  yl 


(4.41) 


H’  =  0 

y1 


(4.42) 


H'  =  — —  sin  0T  E  _ 

el  I  yi 


(4.43) 


4.  The  Perpendicular  Reflected  Wave 

One  of  the  reflected  waves  has  fin  E-vector  polarized  in  the  y' 
direction  E'  .  Equation  (4.8)  in  the  primed  coordinate  system  reads  as 

’  y* 


H  =  -  (sin  0,,0,-cos  0  )  x  (0,  E’  0) 

~x  Hi0  1  1  yx’ 


(4.44) 


and  therefore 


H'  =  E'  cos  0, 
x±  ton0  yx  1 


(4.45) 


H'  .  =  0 

yA 


(4.46) 
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(4.47) 


H*  =  —  sin  0,  E 


z*  CO(i 


i  yi 


5 •  The  Parallel  Reflected  Wave 

In  this  case  the  H-vector  has  only  one  component  H 


Equation  (4.9)  in  the  primed  coordinate  system  reads 


i 

yii 


as 


E  =  - 
~n 


-  (sin  9  0,  -cos  9  )  x  (0,  Hy,|  ,0) 

/.Vi  A  A 


up, 


(4.48) 


and  the  E-field  is  therefore  given  by 


E' 

XII 


cos  0  H' 
i  yii 


(4.49) 


E' 

yn 


=  0 


(4.50) 


E' 


Zll 


— —  sin  0j 
^0 


H' 

yn 


(4.51) 


Now  that  the  field  components  parallel  to  the  boundary  have  been 

evaluated,  they  can  be  substituted  into  the  boundary  condition  equations 

to  obtain  the  magnitude  of  the  generated  waves  as  a  function  of  the 

incident  wave  amplitude.  Continuity  of  the  E',  E',  H',  and  H'  field 

X  y  x  y 
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components  at  the  boundary  give  the  following  four  equations  : 


cos  9  H'  = 

i  y|l 


m  cos  0  /k,  cosY.  sine^  -Uju. 

—  E"  4.  1-5 - s - 5  -  -2-  |h" 

.  ^  \  “Vi  slnVb  VlnV  J 


sin  Y 


E*  +  E' 

yi  yx 


n  sin©  — f,  cos© 


sin  Y 


E"  - 
ya 


mc4if 


kbsin  \ 


H". 

yb 


”  (EyI_  Ey“) 


H' 

y|| 


k  cos0  (n  sin0  -L  cos0  ) 
“  av  a  a' 


m  cose 


sinY 


E"  + 


b  H" 


ya  sin  Y.  yb 


k  n  sine,  -  t  cos0. 

_ 2_  m  E..  _ b _ b  „ 

iopQ  sin^a  ya  MnYb  yb 


Solving  these  equations  yields  the  result 


(4.52) 


E"  =  ~  E* 
ya  yi 


H" 

yb 


=  -  ~  E' 
d.  yi 


(4.53) 


where 


d.  = 


a  b 


{(» 


1  -  slnY.sinV  V"  “‘"V4  cosV  <N.  »°»Vco,ei) 


bsln9b  l 


N, 


-  +  n  sine,  cose 


^cosOj  _  i  cose 


icos0 


(4.54) 


+  m2(cos0  N  +cos6  )  (cose  +N  cose  )> , 
b  d  x  a  a  I  / 


d2  ~  “  -,JT1  Y, 


2cos0,  ,N,  cosY.  sine  , 

I  I  b  b  b  £  .  _  „ 

< - - -  -  —  +  n  sine.  cos0  cose 

(  L,  b  I 


iC098b}  ■ 


(4.55) 
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(cos©a  +  Ng  cosGj)  , 


2  cos©,  m 


d3  =  " 


%  SinYa 


Na  "  kac/u'  and  Nb  “  khc/u 


(4.56) 


E.  BOOKER  QUARTIC 

This  section  is  concerned  with  the  determination  of  N  .  N 

a  b 9 

©a,  and  ©b  appearing  In  Eqs.  (4.53-4.56).  From  Eq .  (4.23),  It  follows 
that 


(4.57) 


Snell's  law  gives  us  ©.  In  the  form 

D 


sin©  slnQ 

sin©  =  - 1  -  _ I 

a  N  v£" 

a  v  1 


(4.58) 


The  corresponding  values  for  the  slow  wave  follow  from  a  simplified 
version  of  the  Booker  quartic  method  (Budden,  1966)  .  The  dispersion 
equation  for  this  wave,  Eq  (4.25),  can  be  rewritten  in  the  form 


Nb  tLl  +  (VLi)cos2*]  -  L0Lx 

Substituting  the  relations 

Nb  cosY  =  l  sin©  +  n  (N.  cos©  ) 
i  d  a 


(4.59) 


(4.60) 


Nb  =  sinGj  +  (Nb  cos©b)‘ 


(4.61) 


into  Eq .  (4.59)  yields  the  equation 


V,1,,2V(Nbc„«V2]  -  »inel+n(Nbco,es)  ]2  =  koKi  (4. 


62) 
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which  is  quadratic  in  N.cosO.  and  has  the  solution 

D  D 


N,cos0,  = 
b  b 


(Lj-Lq)^  sinO^^^-L,^  [n2Lo-sin0I(n2+-t2)+L12(Lo-sin20I)] 

V^W 


(4.63) 


N.  is  then  determined  from  Eq .  (4.61). 
b 

F.  DETERMINATION  OF  FIELD  COMPONENTS  IN  THE  UNPRIMED  COORDINATE  SYSTEM 

In  order  to  carry  out  this  computation  the  components  of  E  ^  and 

H  "  in  the  original  unprimed  coordinate  system  must  be  determined  . 

yb 

This  is  done  by  using  Eq .  (4.3)  to  transform  the  primed  components  of 
the  previous  section  into  the  unprimed  system.  Starting  with  Eqs.  (4.28- 
4.30)  and  (4.35-4.37),  the  following  equations  are  obtained 


E  =  p  E" 
xa  xa  ya 


E  =  p  E"  E  =  p  E" 

ya  xa  ya  za  za  ya 


(4 .64) 


H  .  =  p  .  H". 
xb  xb  yb 


H  ,  =  p  .  H". 
yb  yb  yb 


H  ,  =  P  ,H" 
zb  zb  yb 


(4.65) 


where 


xa 


k  coscpsinO 
_y _ _ _ a 

k  sinY 

a 


ya 


sinf 


L  coscpsinO 
x  T  a 


-  cos9asincp 


za 


kySincpsin©^ 
k  sinY 

a 


xb 


k^coscpsin©^ 

k  sinY 
^  a 


(4.66) 


(4.67) 


(4.68) 


(4.69) 
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[k^costpsin©  1 

■  -  — —  ■  -  cos6  sincp 

kx  *  -I 


kySincpsin©b 

kxain\ 


(4.70) 


(4.71) 


G.  FINAL  SOLUTION  FOR  THE  FAST  ALFVEN  WAVE 

First  it  is  necessary  to  solve  for  the  amplitude  of  the  electric 
field  in  the  x-direction  for  the  fast  wave.  Combining  Eqs.  (4.21), 
(4.22),  (4.53).  and  (4.64),  we  obtain 


r  d2  jkA<4i0M 

E*.<£)  =  -  /<>„,  -  — JT—  d\  <“-72> 

di  8n  k«o 

where  satisfies  Eq.  (4.23). 

If  it  is  assumed  that  (kr)  is  very  large,  it  is  possible  to  apply  the 
method  of  stationary  phase  along  the  lines  given  by  Lighthill  [I960], 
Following  Lighthill f  it  can  be  shown  that  an  integral  of  the  form 

E(r)  =  /F(k)  exp  (-jjk-r)  dk  dk  (4.73) 

✓  ^  y 

when  k  satisfies  the  dispersion  relation 


D(k)  =  0 


can  be  evaluated  in  the  form 


~  F(k>  e*P  (-Jk-r) 
rvfi 


(4.74) 


(4.75* 
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where  7^  is  the  gradient  with  respect  to  k  and  K  is  the  Gaussian  curvature 
of  the  surface 

U(k)  =  0  .  (4.76) 

It  follows  from  Eq.  (4.23)  and  the 
for  K  that 

D  = 

7  D  =  2k 
k  ~ 

!vi  = 2k 

2 

K  =  (4.77) 

u  1^ 


expression  given  by  Lighthill  [1960] 


■>> 


c 


Therefore  the  final  result  for  E  can  be  written  as 

xa 


xa 


‘2  J»2kiU0M^T 


4jtk  „rc 


exp[-j(kzQa+Jk-£)  ]  (4.78) 


where 


k 


(4.79) 


E  and  E  are  found  by  merely  replacing  p  by  p  and  p 
ya  yb  J  xa  ya  Kza  ' 

respectively. 
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H.  FINAL  SOLUTION  FOR  THE  SI/)W  ALFVEN  WAVE 


Proceeding  in  the  same  manner  as  in  the  previous  section,  it  may  be 
shown  that  the  field  for  the  slow  wave  is  given  by 


d3  M 

■Wl>  -  -J\b  “»  dky  <4-> 


where  k  satisfies  Eq.  (4.25).  It  foJ’ows  from  Eq.  (4.25)  that  D  may  be 


written 


D  =  (L^Lj)  k2co,2Y  +  L,k2  -  SL 

c 


=  <Lo_Ll)(Sin  9  kx  +  COS  9  kz)2  +  Llfc2  "  ^  L0L1 

c  1 


(4.81) 


A  straightforward  calculation  then  shows  that 


V  (Lo~Li)  cos  Vsln  9  4  cos  9  +  Lxk/k 

Ivl  [(L0+L1)(L0-L1)cos2  Yb  +  Ll2]1/2 


(4.82) 


Following  Lighthill  [I960],  it  can  also  be  shown  that  the  Gaussian 
curvature  is  given  by 


3  2  2  <  2 

L1  L0  u/c 


k4[LQ2  cos\  +  Lx2  sin\]2 


(4.83) 


It  then  follows  from  Eqs.  (4.75)  and  (4.80)  that,  in  the  asymptotic  limit, 


/  V>\  Pxbd3  Jkik2j;UoM(L02cos2'f,h  +  Sln  . 

(r  fir-,  I  *  '  — - 2 — 2~ W - 1 - ]  • 

V  IVl/  dl  4nk  L, 2' 2L.r 


zO  1  0 


(4.84) 


Again,  H^b  and  Hzb  are  found  by  replacing  Pxb  by  p^b  and  p^b  ,  respectively. 


V.  SOME  PRACTICAL  CONSIDERATIONS 


The  construction  of  a  ground-based  current  loop  large  enough  to 
generate  hydromagnetic  waves  in  the  overhead  ionosphere  involves  certain 
problems  of  a  practical  nature.  What  kind  of  wire  should  be  used  for  the 
loop,  for  example?  Or,  should  the  wire  be  raised  above  the  ground  on 
poles?  Similarly,  the  loop  requires  considerable  power  while  in  operation 
and  the  handling  of  this  power  introduces  further  problems.  Some  of  these 
practical  details  are  discussed  in  this  chapter.  The  object  is  to  show 
that  construction  and  operation  of  the  current  loop  is  feasible  with 
present  technology.  A  loop  with  a  radius  of  10  km  and  carrying  a  current 
of  1000  amps  is  made  the  base  of  our  considerations.  Larger  loops,  with 
greater  current-carrying  capability,  are  feasible  and  will  also  be  dis¬ 
cussed.  However,  it  is  convenient  t.o  have  some  standard  and  the  smaller 
loop  serves  that  purpose. 

A.  DETAILS  OF  THE  LOOP 

1 .  Construction 

It  is  essential  that  the  total  resistance  of  the  current  loop,  and 
therefore  the  resistance  of  the  wire  comprising  it,  be  kept  as  low  as 
possible.  In  other  words,  the  loop  should  be  constructed  of  thick  copper 
or  aluminum  wire.  There  are  two  major  reasons  for  this  requirement. 

First,  the  wire  must  be  capable  of  carrying  Bt  lerst  1000  amp  without 
overheating  and,  second,  the  power  needed  to  drive  a  gi’-en  current  in 
the  loop  is  directly  proportional  to  the  loop  resistance  and  can  only  be 
kept  small  by  minimizing  the  resistance.  However,  wire  resistance  is 
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inversely  related  to  its  diameter,  whereas  its  cost  is  approximately 
directly  related  to  the  diameter.  It  can  be  seen,  therefore,  that 
there  is  a  trade-off  between  the  wire  (and  installation)  costs  and  the 
running  (and  power  conversion)  costs  for  the  current  loop.  We  feel  tnat 
the  benefits  accruing  from  an  increase  in  wire  size  generally  outweigh 
the  added  cost  of  the  wire.  Not  only  are  the  running  and  power  conversion 
costs  reduced  but,  in  addition,  greater  maximum  currents  are  made  possible. 

The  largest  available  wires  are  mostly  constructed  of  aluminum  or 
of  aluminum  over  a  steel  core  (the  aluminum  may  be  either  stranded  or 
unstranded) .  We  would  recommend  stranded  aluminum  over  a  steel  core  of 
gauge  2000  mem  (diameter  w  1-3/4  inches)  and  with  a  resistance  of  approxi¬ 
mately  0 «03  il/km  at  50°C.  The  total  resistance  of  a  10  km  radius  loop 
is  then  just  under  20  and  the  power  required  to  drive  1000  amp  is 

2  Megawatts.  If  the  recommended  wire  is  hung  off  the  ground,  i  .e . ,  in 

free  air,  it  can  carry  a  maximum  current  of  about  3000  amps.  For  this 
large  current  there  would  be  a  rise  of  about  60°C  in  the  temperature  of 
the  wire  above  the  ambient  temperature. 

Insulated  wires  are  available  at  increased  cost  but  have  the  dis¬ 
advantage  in  our  case  that  they  do  not  dissipate  heat  to  the  surrounding  air 
so  readily.  Thus,  their  maximum  current-carrying  capability  is  less  than 
that  of  uninsulated  wire  of  the  same  gauge.  Insulated  wire  has  one 
advantage  in  that  at  a  suitable  location  it  could  conceivably  be  laid  on 
the  ground.  However,  this  further  reduces  the  ability  of  the  wire  to 
dissipate  heat.  We  consider  the  most  effective  loop  to  be  one  constructed 
of  uninsulated  wire  and  strung  on  poles  above  the  ground.  The  poles  add 
extra  expense  but  are  almost  unavoidable.  For  a  loop  of  10  km  radius  or 
larger  there  will  always  be  roads,  streams  and  other  irregularities  of 
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terrain  to  be  crossed,  and  if  the  loop  was  constructed  above  ice  or  snow 
it  would  have  to  be  supported  to  prevent  it  from  sinking.  The  available 
information  suggest?  that  three  turns  of  2000  mcm  gauge  wire  is  probably 
the  maximum  number  of  turns  that  can  be  supported  on  inexpensive  wooden 
poles . 

2 .  Size 

In  principle  there  is  no  limit  on  the  radius  of  the  loop,  and 

there  is  much  to  be  gained  by  making  the  loop  as  large  as  possible. 

Provided  the  loop  can  still  be  treated  as  a  dipole,  its  moment  and 

magnetic  field  vary  as  the  square  of  the  radius,  whereas  the  circumference 

(length  of  wire)  varies  with  the  radius  alone.  Thus,  doubling  the  radius 

of  the  loop  increases  the  field  four  times  (for  a  given  current)  but  only 

requires  as  much  extra  wire  as  in  the  original  loop.  Once  the  dimensions 

start  to  become  comparable  to  the  field  point  distance,  and  the  loop 

can  no  longer  be  treated  as  a  dipole,  increases  in  the  loop  radius  become 

2 

less  effective .  The  change  is  caused  by  the  now  non-negligible  R  term 

in  the  denominator  of  the  equation  for  the  axial  magnetic  field  (see 

equation  A, 3  in  Appendix  A).  Figure  5.1  shows  how  the  change  takes  place. 

Shown  in  the  figure  is  the  variation  of  the  axial  magnetic  field  at  100  km 

11  2 

for  a  loop  of  constant  magnetic  momen t  (it.  10  amp.  m  )  but  with  a  radius 
varying  from  0  km  (i.e.  dipole)  to  150  km.  It  can  be  seen  that  a  rapid 
decrease  in  magnetic  field  commences  when  the  radius  reaches  about  20  km. 
The  rate  of  decrease  is  greatest  in  the  radius  range  20  -  50  km  and 
becomes  less  beyond  the  50  km  point  of  inflection.  However,  it  should  be 
pointed  out  that  the  field  produced  by  the  50  km  loop  ie  only  about  one 
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third  less  than  the  dipole  field,  whereas  the  field  of  the  100  km  loop 
is  about  half  the  field  of  the  50  km  loop.  In  practice,  of  course,  the 
magnetic  moment  of  the  loop  increases  with  increasing  radius  and  greater 
axial  magnetic  fields  can  still  be  produced  by  keeping  the  current  unchanged 
and  increasing  the  radius  of  the  loop.  The  point  is  that  doubling  the  loop 
radius  no  longer  necessarily  increases  the  axial  field  four  times.  For 
example,  given  a  certain  fixed  loop  current,  an  increase  of  the  loop  radius 
from  50  km  to  100  km  only  doubles  the  axial  magnetic  field  at  100  km. 

With  the  above  considerations  in  mind,  we  would  limit  the  radius  of 

the  loop  to  a  value  in  the  range  25  -  50  km.  The  field  produced  by  such 

a  loop  can  be  quite  substantial  at  E- region  heights.  A  three  turn  loop 

of  radius  50  km  carrying  a  current  of  3000  amps  has  a  resistance  of  about 

28ft  (2000  mem  aluminum  wire),  a  loop  moment  of  225jt  x  1011  amp.m2 
13  2 

(7.1  x  10  fimp .  m  )  and,  at  an  altitude  of  100  km  above  the  loop,  produces 
an  axial  field  of  close  to  one  gamma.  This  amplitude  is  large  compared 
with  typical  magnetic  field  amplitudes  usually  associated  with  micropulsations. 
The  field  can  be  increased  still  further  but  with  somewhat  greater  effort 
and  less  efficiency.  We  consider  the  three-turn  loop  of  50  km  radius,  with 
a  maximum  current  of  3000  amps,  to  be  the  limit  to  what  might  be  comparatively 
easily  achieved.  (The  power  requirement  for  such  a  loop  is  huge,  about 
250  Mw,  but  huge  powers  are  available  from  the  utilities  and  we  are  here 
investigating  the  limits  to  a  practical  system)  . 

In  the  remainder  of  this  chapter  the  three-turn  50  km  radius  loop, 
carrying  3000  amp,  will  be  referred  to  as  "the  50  km  loop”.  The  standard 
single-turn,  10  km  radius  loop,  carrying  1000  amp,  will  similarly  be 
referred  to  as  ’’the  10  km  loop” . 
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LOOP  RADIUS  (Km) 


Variation  of  the  self-inductance  of  a  single 
circular  loop,  constructed  of  wire  of  gauge 
2000  racm  (diameter  1.75  inches),  as  the  loop 
radius  is  increased  from  about  0  km  to  100  km. 


3.  Self-Inductance 


The  self-inductance  of  a  single  circular  loop  of  radius  R  and 
wire  radius  r  is  given  by 


L  =  hqR  [ln(2|>  -  I]  .  (5.1) 

There  are  some  approximations  involved  in  deriving  this  formula,  but 
they  all  appear  to  be  satisfied  in  this  case  (r  «  R,  as  required,  and 
the  skin  effect  will  be  negligible  at  1  Hz) .  Figure  5  .2  shows  the 
variation  of  L  for  loops  varying  in  radius  from  ~  0  km  to  100  km  and 
constructed  of  wire  of  gauge  2000  mcm.  It  can  be  seen  that  L  is  roughly 
proportional  to  R  and  varies  from  168  mH  for  a  10  km  loop  to  1960  mH 
for  a  100  km  loop.  The  inductance  values  are  not  exceptionally  large. 

To  illustrate,  the  multi-turn  solenoids  used  to  record  micropulsations  at 
Stanford  are  about  5  feet  long  and  6  inches  in  diameter  and  have  inductances 
of  ~  16  H. 

The  importance  of  the  above  self-inductances  can  be  assessed  by 
comparing  the  ohmic  resistance  of  the  coils  with  their  inductive  reactances 
uL  .  With  2000  mcm  aluminum  wire,  the  resistance  of  the  loop  varies  from 
about  2Q  for  the  10  km  loop  to  9 .4fi  for  a  50  km  single-turn  loop.  The 
corresponding  inductive  reactances  at  1  Hz  vary  from  about  1  fi  to 
about  G  n  and  are  always  less  than  the  ohmic  resistances.  These  data 
indicate  that  the  self- inductance  of  the  loops  will  not  interfere  signifi¬ 
cantly  with  their  use  as  micropulsation  generators.  On  the  other  hand,  the 
self-inductance  effects  will  not  be  completely  negligible  and  should  be 
carefully  evaluated  in  any  design  program. 
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Figure  5.3  Block  diagram  of  a  power  system 
for  the  10  km  loop. 
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B.  POWER  MANAGEMENT  AND  FREQUENCY  GENERATION 

Early  In  this  study  it  became  clear  that  there  were  several 
alternative  ways  to  produce  a  large  alternating  current  in  a  conducting 
loop.  The  method  we  chose  to  investigate  in  some  detail  is  probably  the 
most  straightforward  of  these  methods:  the  current  is  either  reversed  or 
switched  on  and  off  repetitively  by  suitable  control  devices.  A  "square 
wave"  magnetic  field  variation  is  produced  with  a  frequency  that  is 
dependent  upon  the  switching  rate.  We  do  not  assert  that  this  method  is 
necessarily  the  best  for  the  production  of  the  required  large  alternating 
current  in  the  loop.  At  the  same  time,  the  switching  method  does  have 
a  number  of  desirable  features  and  is  worth  serious  consideration  for 
artificial  micropulsation  generation. 

The  major  disadvantage  of  a  "square  wave"  magnetic  field  variation 
is  the  presence  of  harmonics.  These  harmonics  absorb  roughly  one  quarter 
(23.4%)  of  the  input  power  and  the  higher  frequencies  (~  20  Hz)  may  also 
have  undesirable  interference  effects.  The  overall  loss  of  power  is 
probably  the  most  severe  problem.  Comparatively  little  power  is  absorbed 
by  the  higher  frequency  harmonics,  compared  to  the  lower  frequency 
harmonics,  and  they  are  more  severely  attenuated  by  the  earth  and  ionosphere. 
Their  possible  interference  effects  are  thus  restricted  to  an  area  close 
to  the  current  loop. 

Figure  5.3  shows  a  block  diagram  of  the  power  system  for  the 
standard  10  km  loop.  A  voltage  of  close  to  2  kV  is  required  to  drive 
the  maximum  loop  current  of  1000  amps  and  it  is  obtained  with  little 
difficulty  from  the  standard  60  -  70  kV  power  transmission  voltage.  As 
shown  in  the  block  diagram  the  voltage  entering  the  system  is  first 
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transformed  clown  to  a  lower  value  (~  4kV)  in  a  switchyard  section.  The 
switchyard  also  serves  to  protect  the  power  company  and  the  following 
converter  stage.  In  the  converter  stage  the  ~  4kV  three  phase  A ,C . 
voltage  is  converted  to  the  desired  switched  D.C.  voltage  that  is 
applied  directly  to  the  current  loop  "antenna".  Transformers  in  the 
converter  section  adjust  the  A.C.  voltage  to  give  the  required  2kV  D.C. 
voltage  amplitude.  The  control  circuitry  necessary  for  all  the  operations 
in  the  converter  stage  is  indicated  by  a  separate  block  in  the  diagram. 

The  resistance  of  the  50  km  loop  is  about  28  H  and  a  voltage  of 
close  to  80  kV  is  required  to  drive  the  desired  maximum  current  of 
3000  amps.  It  would  obviously  be  very  difficult  to  supply  all  the 
voltage  and  power  at  one  point  on  the  loop.  A  possible  solution  is  to 
space  several  substations  around  the  circumference  of  the  loop  with 
each  substation  supplying  part  of  the  required  voltage  and  having  a 
power  system  similar  to  the  one  shown  in  Figure  5.3. 

C .  LOCATION  OF  THE  CURRENT  LOOP 

The  most  important  positive  factor  in  the  choice  of  a  site  for  the 
current  loop,  and  we  will  initiallv  restrict  consideration  to  sites  within 
the  continental  United  States,  is  either  a  surface  layer  of  low  to  medium 
electrical  conductivity  or  at  least  onli  a  thin  surface  layer  of  high 
conductivity.  Other  important  factors  are  (a)  suitable  location  relative 
to  the  geomagnetic  field;  (b)  remoteness  of  the  site  (to  minimize  inter¬ 
ference  with  telephone  and  power  circuits);  (c)  power  availability; 

(d)  security  (to  prevent  theft  of  wire  and  other  components);  and  (e)  low 
cost . 
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Naturally-occurring  Pc  1  micropulsations  usually  originate  on 
geomagnetic  field  lines  that  extend  out  to  about  4  earth  radii  or  more, 
at  their  equatorial  points,  from  the  earth's  center.  For  the  current 
loop  to  be  suitably  located  relative  to  these  field  lines  it  would  have 
to  be  constructed  in  Canada  or  Alaska.  While  we  have  no  specific  objection 
to  either  location,  it  is  clear  that  each  has  a  number  of  minor  disadvantages 
and  no  major  advantages  other  than  the  one  discussed.  Furthermore,  it  can 
be  argued  that  an  artificial  source  of  micropulsations  does  not  necessarily 
have  to  be  located  in  the  most  favorable  location  for  the  generation  of 
the  natural  signals.  On  balance,  therefore,  we  decided  that  a  site  south 
of  the  Canadian- USA  border  would  be  adequate  for  the  current  loop. 

Factors  (b)  -  (e)  all  suggest  that  the  site  should  be  in  the 
Western  part  of  the  United  States  where  there  are  large  expanses  of 
largely  uninhabited  desert  (often  in  restricted  areas)  and  where  there 
are  large  hydroelectric  power  stations  capable  of  supplying  abundant 
power.  However,  the  surface  conductivity  of  the  desert  areas  is  never 
very  low  and,  in  some  locations,  is  undesirably  high.  Keller  et  al.  (1966) 
show  a  surface  resistivity  map  of  the  United  States  (Figure  5.4),  and  the 
three  areas  of  lowest  conductivity  are  in  the  New  England  -  Appalachia 
region,  in  northern  Minnesota,  and  in  a  section  of  the  northern  Rocky 
Mountains.  In  the  west  and  mid-west  the  conductivity  is  predominantly 
quite  high.  The  areas  of  moderate  conductivity  in  this  region  occur 
irregularly  between  the  Sierra  Nevada  in  California  (and  the  Cascade 
Range  to  its  north)  and  the  eastern  termination  of  the  Rocky  Mountains. 

Keller  et  al,  do  not  indicate  the  depth  of  the  surface  layers  whose 
resistivities  are  shown  in  the  figure.  For  the  high  conductivity  areas 
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Surface  resistivity  of  the  United  States 
(After  Keller  et  al.,  1966) 


66 


in  the  west  the  depth  probably  varies  between  10  m  and  1  km  [e.g., 
Jackson  (1966)]  , 

The  situation  in  central  Nevada  is  of  ^special  interest,  because 
our  preliminary  investigations  indicate  that  an  A EC  test  site  may 
possibly  be  available  for  construction  of  the  current  loop.  The  site 
in  question  is  restricted  to  authorized  personnel,  is  presently  inactive 
and  has  10  Mw  of  power  capacity  already  available.  The  surface  layer 
conductivity  in  Nevada  varies  from  moderate  to  high,  with  the  high 
conductivities  mostly  confined  to  the  valleys.  Jackson's  (1966)  results 
indicate  that  the  high  conductivity  layers  are  thin  relative  to  their 
skin  depth  at  micropulsation  frequencies.  Thus,  central  Nevada  may 
well  be  a  very  suitable  location  for  construction  of  the  current  loop. 

The  only  locations  outside  the  United  States  that  may  have 
unique  advantages  as  sites  for  the  current  loop  are  Greenland  and 
Antarctica.  Both  continents  have  extremely  thick  ice  sheets  resting 
on  rock.  From  the  data  in  Chapter  3,  the  ice  should  have  very  low 
conductivity,  and  if  the  underlying  rock  is  also  of  low  conductivity 
(i.e.,  non-sedimentary  rock)  either  location  could  be  suitable  for  the 
current  loop.  Other  advantages  are  (a)  excellent  coverage  of  desirable 
field  lines,  (b)  practically  unlimited  level  area  for  construction 
of  the  loop,  with  no  right-of-way  problems  (arrangements  would  have  to 
be  made  with  the  appropriate  authorities  for  the  use  of  either  site), 

(c)  few  people,  and  (d)  little  opportunity  for  interference  with  power 
and  telephone  circuits.  However,  while  all  these  advantages  are  highly 
desirable,  they  are  not  crucial  for  the  success  of  the  current  loop 
experiment.  Greenland  and  Antarctica  have  one  very  great  disadvantage 
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at  the  present  time:  inadequate  power.  Until  more  power  can  be 
provided,  or  the  current  loop  can  be  made  to  operate  on  powers  at  the 
kilowatt  level  (i.e.,  non-continuous  operation),  we  do  not  consider 
Greenland  or  Antarctica  to  be  feasible  sites  for  the  current  loop. 

D.  COST  OF  A  GROUND-BASED  CURRENT  LOOP  SYSTEM 

We  have  obtained  preliminary  estimates  from  several  companies 
experienced  in  heavy  electrical  engineering  for  construction  of  the 
10  km  loop  and  its  power  system.  Very  briefly,  the  cost  of  the  loop 
itself  installed  on  poles  varies  from  $360,000  to  ~  $500,000.  Control 
cabling,  lightning  cable,  and  three-phase  A.C.  lines  cost  ~  $100,000 
and  the  D.C.  converter  stage  is  ~  $200,000.  (Two  converters  are  needed, 
one  for  the  positive  and  another  for  the  negative  voltages,  and  their 
cost  is  ~  $50/Kw) .  Allowing  for  labor  and  contingencies,  the  cost  of 
the  10  km  loop  and  power  system  would  probably  be  close  to  $1  million. 

Projecting  the  estimates  received  for  the  10  km  loop,  we  estimate 
a  cost  of  ~  $12  million  for  the  50  km  loop  and  power  system. 

These  estimates  are  all  preliminary  and  subject  to  change.  They  are 
included  for  completeness  and  to  give  perspective  to  the  concepts  and 
calculations  in  this  report. 

E.  ENVIRONMENTAL  ASPECTS  OF  LARGE  CURRENT  LOOPS 

It  is  probably  unnecessary  to  point  out  that  current  loop  of  the 
size  considered  in  this  report  will  only  be  possible  in  areas  that  are 
largely  uninhabited.  Nevertheless,  particular  concern  must  be  paid  to 
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to  the  effect  of  the  magnetic  field  on  both  human-  and  wild- life,  on 
plants,  and  on  telephone  and  power-line  circuits.  Fortunately,  the 
magnetic  field  from  a  current  loop  has  several  features  that  limit  its 
environmental  impact.  Perhaps  the  most  important  of  these  features  is 
a  very  rapid  drop-off  in  the  amplitude  of  the  field  with  distance.  In 
addition,  it  is  advantageous  that  the  field  amplitudes  are  mostly  small 
relative  to  the  geomagnetic  field  and  that  their  frequency  of  oscillation 
is  very  low . 

At  distances  large  compared  to  the  radius  of  the  loop  the  dipole 
approximation  becomes  valid  and  the  field  decreases  as  the  inverse 
cube  of  distance.  Calculations  with  the  finite  loop  expressions  given 
by  Smythe  (Appendix  A;  see  Figure  5.5)  indicate  that  in  the  plane  of  the 
loop  the  dipole  approximation  is  accurate  to  within  about  5  percent  at 
a  distance  of  four  times  the  radius  from  the  center  of  the  loop.  (Note 
that  the  plane  of  the  loop  represents  the  earth's  surface).  For  the 
50  km  loop,  the  field  at  200  km  (i.e.,  at  4R)  from  the  center  of  the  loop 
is  close  to  ly .  Magnetic  variations  of  this  amplitude  and  with  micro¬ 
pulsation  frequencies  are  not  uncommon  in  the  auroral  zones  and  appear 
to  have  no  harmful  effects  whatsoever.  We  will  refer  to  the  distance 
for  a  ly  field  as  the  "critical  distance". 

For  distances  greater  than  the  critical  distance  the  current  loop 
is  likely  to  have  a  negligible  effect  on  everything  except  sensitive 
magnetometers.  Our  calculations  indicate  that  the  critical  distance  will 
vary  from  approximately  200  km  (50  km  loop)  to  about  30  km  (10  km  loop) 
for  the  loop  sizes  that  might  be  used  in  micropulsation  generation. 
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MAGNETIC  FIELD  IN  THE  PLANE  OF  THE  LOOP  (gamma) 


| 

j 


DISTANCE  FROM  CENTER  (km) 


Figure  5.5  Magnetic  field  variations  in  the  planes  of  the  two 

current  loops  discussed  in  the  text.  The  field  inside 
et’ch  of  the  loops  is  opposite  to  the  field  outside  and 
its  variation  is  indicated  by  a  dashed  line.  For 
comparison,  the  equivalent  dipole  field  variations  are 
also  shown. 
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Figure  5.5  compares  the  magnetic  fields  produced  in  the  z  =  0  plane 
(i.e.,  in  the  plane  of  each  loop)  by  the  10  km  and  50  km  loops.  The 
equivalent  dipole  fields  are  shown  together  with  the  finite  loop  fields 
calculated  by  using  equations  (A. 5)  -  (A. 7)  in  Appendix  A.  The  magnetic 
field  from  a  current  loop  is  perpendicular  to  the  plane  of  the  loop  at 
all  points,  but  there  is  a  180°  change  of  direction  if  the  field  at 
points  inside  and  outside  the  loop  are  compared.  Figure  5.5  shows  that 
the  largest  fields  produced  by  the  loop  occur  at  the  loop  surfaces.  They 
do  not  rise  to  infinity,  as  do  the  dipole  fields  at  the  origin,  but  reach 
values  that  may  be  obtained  from  Equation  (A. 8)  in  Appendix  A.  Except  in 
the  immediate  vicinity  of  the  loops  themselves  (within  ~  12  m  for  the 
50  km  loop)  the  field  amplitudes  are  all  less  than  the  earth's  magnetic 
field.  It  is  interesting  to  compare  the  fields  at  the  center  of  the  two 
loops.  The  field  at  the  center  of  the  10  km  loop  is  62 .8y,  or  about  half 
the  113. ly  field  at  the  center  of  the  50  km  loop.  There  is  a  surprisingly 
small  difference  in  these  fields  considering  that  the  moment  of  the  50  km 
loop  is  225  times  as  large  as  the  moment  of  the  10  km  loop.  However,  the 
field  at  the  center  of  a  current  loop  depends  on  i/R  (see  Equation  (A .4)  in 
Appendix  A)  and  it  can  be  seen  that  a  simultaneous  increase  of  both  the 
loop  current  and  radius  does  not  necessarily  increase  the  central  field, 
even  though  the  moment  of  the  loop  and  its  field  at  large  distances  may 
be  increased  substantially. 

For  distances  less  than  the  critical  distance  it  is  possible  that 
significant  voltages  will  be  generated  in  telephone  or  power  circuits  by 
the  alternating  magnetic  field  of  a  current  loop.  A  necessary  condition, 
of  course,  is  that  the  circuits  should  form  closed  loops  with  sufficient 
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area  inside  the  circle  defined  by  the  critical  distance.  The  limited 
area  of  strong  magnetic  field,  the  low  frequencies  and  the  location  of 
the  current  loop  within  an  uninhabited  area  will  greatly  reduce  the 
possibility  r"c  interference  with  telephone  and  power  circuits  .  Telephone 
and  power  lines  to  facilities  close  to  the  current  loop,  and  terminating 
at  those  facilities,  will  very  probably  be  unaffected. 

Very  little  information  is  available  on  the  effect  of  weak  to 
moderate  alternating  magnetic  fields  (amplitudes  ^  100  gauss)  on  plant 
and  animal  life.  However,  there  are  two  good  reasons  to  believe  that  the 
magnetic  fields  produced  by  the  current  loop  discussed  in  this  report 
will  have  no  harmful  effects.  First,  the  amplitudes  of  the  magnetic 
field  are  nearly  everywhere  very  much  smaller  than  the  earth's  magnetic 
field  .  They  are  many  orders  of  magnitude  smaller  than  some  of  the 
magnetic  fields  commonly  encountered  in  our  society,  and  repeated  exposure 
of  human  life  to  these  strong  fields  appears  to  have  negligible  effect. 
Second,  the  comparatively  strong  fields  existing  close  to  the  loop  are 
not  unique  to  these  systems.  Magnetic  fields  of  similar  strength  must 
surround  the  larger  utility  lines  and  yet  these  fields  appear  to  be  quite 
harmless.  In  addition,  the  frequency  of  variation  of  these  utility  fields  is 
considerably  higher  (~  60  times)  than  the  frequencies  under  consideration 
for  the  current  loop.  Assuming  that  the  effect  of  alternating  magnetic 
fields  on  plant  and  animal  life  is  proportional  to  the  induced  currents 
generated  by  the  fields,  the  effects  created  by  the  current  loops,  if  any, 
must  be  about  sixty  times  less  than  those  produced  by  the  larger  utility 
lines.  Obviously,  further  consideration  should  be  given  to  the  biological 
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effects  of  weak  to  moderate  magnetic  fields  if  a  large  current  loop  is 
ever  to  be  constructed,  but  at  this  stage  we  feel  that  the  magnetic  field 
of  such  a  loop  is  entirely  innocuous . 


VI .  CONCLUSION 


The  results  presented  in  the  previous  chapters  indicate  that  it  is 
technically  feasible  to  construct  a  current  loop  on  the  ground  with  a 
magnetic  moment  as  large  as  7  x  10  amp.m  and  that  such  a  loop  will 
produce  a  steady  magnetic  field  of  about  1  gamma  in  the  E-region  of  the 
ionosphere.  Furthermore,  this  field  can  be  alternated  at  micropulsation 
frequencies  by  switching  the  loop  current.  Some  of  the  power  input  to  the 
loop  is  dissipated  through  the  generation  of  harmonics,  and  earth  currents 
are  induced  that  reduce  the  magnetic  field  of  the  loop.  However,  at 
micropulsation  frequencies,  it  still  appears  possible  to  produce  a  magnetic 
field  variation  in  the  E-region  with  an  amplitude  of  about  l/2  gamma.  An 
artificial  1  Hz  field  variation  of  this  amplitude  is  large  compared  to 
most  natural  Pc  1  micropulsation  amplitudes  observed  beneath  the  ionosphere. 
It  seems  reasonable  to  expect,  on  this  basis,  that  the  alternating  magnetic 
field  from  the  current  loop  could  excite  significant  Pc  1  hydromagnetic 
wave  amplitudes  in  the  lower  ionosphere.  Hm  waves  of  other  micropulsation 
frequencies  might  also  be  excited  by  suitable  changes  in  the  switching 
frequency  of  the  current  in  the  loop. 

The  actual  amplitudes  of  the  two  hydromagnetic  wave  modes  produced 
in  the  E-region  by  the  current  loop,  i.e.,  the  fast  and  slow  Alfven  waves, 
can  be  calculated  by  using  the  appropriate  expressions  derived  in  Chapter  IV. 
(The  magnetic  moment  M  of  the  loop  may  be  replaced  by  an  "effective"  magnetic 
moment  to  allow  for  the  generation  of  harmonics  and  earth  currents).  Thus, 
it  is  possible  to  answer  the  following  question:  can  a  current  loop 
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with  an  alternating  magnetic  moment  of  amplitude  7  x  10  amp.m 
produce  artificial  micropulsations  at  large  distances  from  the  loop? 

To  answer  this  question  completely  the  amplitudes  of  the  two  hm  modes 
produced  in  the  E-region  must  be  calculated  and  their  subsequent  propagation, 
and  propagation  losses,  in  the  ionosphere  and  magnetosphere  considered  in 
some  detail.  The  calculations  are  lengthy  and  we  restricted  ourselves  to 
preliminary  estimates  of  the  amplitudes  of  the  hm  modes.  For  the  hm  mode 
that  is  ducted  in  the  ionosphere  at  Pc  1  frequencies  (i.e.,  the  fast  Alfven 
mode)  our  estimates  indicate  that  mode  amplitudes  of  at  least  1  -  10my 

13 

are  achievable  with  a  ground-based  current  loop  of  maximum  moment  7  x  10 
2 

amp.m  .  Thus,  under  nighttime  conditions,  with  little  absorption,  it 

should  be  possible  to  produce  artificial  Pc  1  micropulsations  with  amplitudes 

greater  than  lmy  at  large  distances  from  the  current  loop.  This  conclusion 

is  independently  verified  by  the  work  of  Greifinger  (1972)  .  According  to 

13  2 

Greifinger,  a  horizontal  magnetic  loop  with  a  moment  of  jt  x  10  amp.m  , 
located  on  a  single- la ye red  earth  of  conductivity  1-2  x  10  4  mho/m,  can 
produce  an  artificial  1  Hz  micropulsation  of  amplitude  ~  7my  at  a  distance 
of  5000  km  from  the  loop.  This  result  applies  only  to  the  hm  mode  that 
is  ducted  in  the  F-region  oi  the  ionosphere  and  for  high-latitude  nighttime 
conditions.  However,  it  confirms  that  artificial  Pc  1  micropulsations 
can  be  generated  by  a  ground-based  current  loop  at  large  distances  from  the 
loop . 

The  theory  and  calculations  included  in  this  report,  and  in  Greifinger 's 
(1972)  work,  need  to  be  extended  in  several  areas  in  order  to  obtain  an 
overall  picture  of  the  capability  of  the  current  loop.  Most  pressing, 
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perhaps,  is  a  complete  calculation  of  the  hm  energy  propagating  along  a 
field  line  to  a  conjugate  point.  We  are  not  quite  as  pessimistic  as 
Greifinger  concerning  the  possibility  of  completing  the  required  calcula¬ 
tion  with  present  knowledge,  but  it  is  undoubtedly  a  complex  and  difficult 
problem  (c.f.  Kitamura  and  Jacobs,  1968).  The  calculations  also  need  to 
be  extended  to  micropulsation  frequencies  other  than  Pc  1.  For  the 
lowest  micropulsation  frequencies  the  single- layered  earth  model  is 
inadequate  and  the  effect  of  the  earth  currents  on  the  performance  of 
the  loop  should  be  evaluated  by  using  2-  and  3-layer  models. 

The  loop  power  system  discussed  in  Chapter  V  is  designed  for  CW 
(continuous-wave)  operation  and  there  are  many  advantages  to  this  method 
of  operation.  For  example,  the  artificial  micropulsations  will  probauly 
be  quite  weak  at  large  distances  from  the  loop  and  a  CW  signal  enables 
integration  to  be  carried  out  in  a  simple  and  direct  manner.  Under 
some  circumstances,  however,  pulsed  operation  of  the  loop  may  have  some 
advantages,  particularly  when  there  is  a  limit  on  the  power  available 
(see  Chapter  V,  Section  C) .  Thus,  further  engineering  studies  on  the  loop 
power  system  would  be  worthwhile. 

There  have  been  two  other  methods  proposed  for  the  generation  of 
artificial  micropulsations:  (a)  by  using  large  loop  antennas  on  satellites 

in  the  ionosphere  and  (b)  by  using  large  VLF  transmitters.  The  first  of 

-3 

these  methods  avoids  the  r  drop-off  in  magnetic  field  that  occurs  between 
a  ground-based  current  loop  and  the  ionosphere.  However,  the  required 
loop  moment  for  the  generation  of  hm  waves  is  probably  still  quite  large 
and  it  is  very  difficult  to  arrange  these  large  moments  for  satellite- 
borne  loops.  The  second  method,  proposed  by  Bell  (1972),  is  based  upon 
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the  precipitation  of  energetic  electrons  from  the  magnetosphere  by  VLF 
wave  injection.  Periodic  transmissions  from  a  VLF  ground  or  satellite 
transmitter  are  used  to  precipitate  periodic  bursts  of  energetic  electrons 
into  the  ionosphere.  These  bursts  induce  periodic  changes  in  the  D-  and 
E-reglon  conductivities,  resulting  in  periodic  fluctuations  in  the  currents 
that  flow  at  these  levels  and  hence  micropulsation-like  magnetic  field 
changes  on  the  ground.  [The  precipitating  electrons  do  not  appear  to 
generate  hm  waves  in  the  ionosphere  and  thus  the  magnetic  field  changes  do 
not  fit  the  definition  adopted  in  this  report  for  artificial  micropulsa¬ 
tions]  .  There  are  some  disadvantages  associated  with  this  method.  The 
periods  are  likely  to  be  restricted  to  values  greater  than  about  15 
seconds,  and  the  ability  of  the  experimenter  to  precipitate  electrons  will 
depend  very  greatly  upon  the  variable  state  of  the  magnetosphere.  Never¬ 
theless,  the  VLF  injection  method  is  interesting. 

Another  possible  method  of  generating  artificial  micropulsations 

is  to  use  an  electric  dipole  instead  of  the  magnetic  dipole  described  in 

this  report.  If  the  moment  of  the  electric  dipole  is  large  enough  and  can 

be  varied  at  micropulsation  frequencies  we  would  expect  hm  waves  to  be 

generated  in  the  overhead  ionosphere.  We  have  not  calculated  the  amplitude 

of  the  hm  modes  generated  in  the  ionosphere  for  an  electric  dipole;  neither 

» 

have  we  considered  the  effect  of  the  earth  on  the  dipole  field.  However, 
we  consider  the  electric  dipole  method  for  generating  micropulsations  less 
likely  to  be  successful  than  the  magnetic  dipole  method  because  the  technology 
for  producing  large  currents  is  much  better  developed  than  the  technology 
for  large  charges  .  An  electric  dipole  is  also  more  susceptible  to  the 
weather  than  a  magnetic  dipole.  Interestingly,  it  has  been  suggested 
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(Fraser-Smith  &  Roxburgh,  1969)  that  the  electromagnetic  fields  produced 
during  lightning  discharges  may  trigger  certain  Pc  1  micropulsations,  and 
it  is  well-known  that  lightning  discharges  are  associated  with  large 
electric  dipole  fields  in  clouds  or  between  clouds  and  the  earth. 

To  sum  up,  it  appears  possible  to  generate  artificial  micropulsatrlons 
with  a  large  ground-based  current  loop.  The  micropulsation  frequencies, 
and  to  some  extent  their  amplitude,  are  under  the  control  of  the  operator 
of  the  loop.  There  are  no  limits  on  the  time  of  day,  or  the  state  of 
geomagnetic  activity,  during  which  the  loop  may  be  used  (although  the 
class  of  micropulsations  that  the  experimenter  might  wish  to  generate 
will  vary  with  time  of  day  and  geomagnetic  activity)  .  We  believe  such 
a  loop  would  significantly  advance  knowledge  of  the  ULF  behavior 
of  the  ionosphere  and  magnetosphere.  Also,  because  there  is  a  limit  to 
what  can  be  achieved  by  passive  observation,  the  future  of  micropulsation 
research  is  clearly  circumscribed  without  a  method  of  controlled  artificial 
generation.  The  construction  of  a  current  loop  generation  system,  and  the 
possibility  of  controlled  experiments,  would  greatly  expand  our  knowledge 
of  micropulsation  phenomena  and  provide  a  powerful  new  tool  for  probing 
both  the  magnetosphere  and  the  interior  of  the  earth. 
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Appendix  A 

MAGNETIC  FIELD  OF  A  CURRENT  LOOP  IN  FREE  SPACE 


1.  Dipole  Expressions.  The  magnetic  field  components  for  a  magnetic 
dipole  of  moment  M  directed  along  the  positive  z  direction  of  a  cylindrical 
system  of  coordinates  are  given  by: 
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where  p  =  (z2  +  r2)  '  is  the  distance  from  the  dipole  to  the  field 


point . 


When  the  field  point  is  on  the  z-axis: 
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2.  Finite-Loop  Expressions .  The  general  expressions  for  the  magnetic 
field  produced  by  a  current  loop  of  finite  size  are  derived  by  Smythe  (1950) . 
They  involve  elliptic  integrals  and  are  quite  complicated  in  their  general 
form.  The  expressions  tc«e  a  simple  form  for  field  points  on  the  axis, 
however,  and  it  is  interesting  to  compare  this  simple  form  with  the 
equivalent  dipole  expressions. 
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For  a  current  loop  of  radius  E  lying  In  the  z  =  0  plane  of  a 
cylindrical  system  of  coordinates  and  with  its  center  at  the  origin  we 
have,  for  points  on  the  z  axis: 
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where  M  =  I  •  «R  is  the  loop  moment  and  it  is  assumed  that  the  vector 
moment  is  directed  along  the  positive  z  direction.  At  the  center  of  the 
loop  we  have 


(A  .4) 


When  the  field  points  are  in  the  plane  of  the  loop  (z  =  0)  the 
general  expressions  given  by  Smythe  (1950)  for  the  magnetic  field  of  the 
loop  reduce  to 
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where  K  and  E  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind,  respectively.  That  is 


with 


K  = 


n/2  -1/2 

f  [1  -  k  sinZ0]  d@ 

rr/2  2  9 

E  =  *'*[1  -  k  sinZ0]  d0  (A. 6) 
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For  r  R,  i.e.,  for  points  very  close  to  the  wire  comprising  the  loop, 
equations  (A. 5)  reduce  to 


B  =  0  =  B - 

r  0 

B  -  ^ 

z  2jt(R-  r) 


(A. 8) 


where  the  equation  for  B  is  now  identical  to  the  equation  for  the 

z 

field  produced  at  distance  Ar  =  (R-r)  by  a  straight  wire  carrying  a 
current  I . 
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3.  Comparison  of  the  Dipole  and  Finite-Loop  Fields .  For  large  values 


of  z,  i.e.,  large  compared  to  the  loop  radius  R,  the  axial  field  of  the 

finite  loop  will  obviously  approximate  the  axial  field  for  a  dipole.  Even 

when  z  is  only  a  little  larger  than  the  loop  radius  it  is  often  still  a 

good  approximation  to  use  the  dipole  expression  for  the  axial  field. 

Consider  a  particular  example  that  is  relevant  to  the  subject  of  this 

report .  Suppose  the  moment  of  both  the  loop  and  the  dipole  is  n  • 1011 
2 

amp.  m  ,  with  the  loop  having  a  radius  of  10  km  and  carrying  a  current 
of  1000  amp.  At  z  =  100  km  (z  =  10R)  the  dipole  field  is  62-8  x  1012wb/m2, 
or  62*8  my.  The  equivalent  field  for  the  loop  is  61 ’9  my.  Thus,  there 
is  little  difference  between  the  two  field  magnitudes  at  a  distance 
equivalent  to  the  height  of  the  ionospheric  E- region. 

The  difference  between  the  two  fields  for  a  loop  of  radius  100  km 
may  also  be  of  interest.  The  loop  in  this  case  can  only  have  a  current 
of  10  amp  if  the  moment  is  to  be  the  same  as  before.  For  this  example 
the  loop  field  at  100  km  (z  =  R)  is  22-2  my,  or  about  one  third  the 
dipole  field. 


APPENDIX  B 


ALTERNATIVE  EXPRESSIONS  FOR  THE  TOTAL  MAGNETIC 
field  PRODUCED  BY  A  GROUND-BASED  CURRENT  LOOP 


It  follows  from  equations  (3.38)  and  (3.45)  that  the  total 
magnetic  field  for  z  0  has  components 
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Consider  the  equation  for  a  in  more  detail.  This  equation  may  be 
written 
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where  equation  (?.0)  for  the  skir  depth  6  has  been  used. 
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The  components  of  total  magnetic  field  become 
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On  the  axis  of  the  loop,  equations  (B.5)  give 
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where 


Equations  (B.6)  show  that  on  the  axis  varies  with  frequency  and 

conductivity  as  (cjcj)  ^  ,  provided  the  approximation  ~ —  »  1 

11  2  6  X 

is  valid.  For  z  =  100  km  and  =  n.lO  amp. m  we  obtain 
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Equations  (B.5)  and  (B.6)  are  convenient  alternative  expressions  for 
the  total  magnetic  field  produced  by  the  ground-based  loop  (with  the 
earth  assumed  to  be  a  single  layer  of  conductivity  a) , 

We  now  investigate  the  conditions  for  equations  (B.5)  and  (B.6) 
to  be  valid.  The  condition  — »  1  may  be  written  \«6-1  .  Suppose 

rr 

is  a  value  of  X  that  can  be  considered  much  smaller  than  6_1,  i  .e . , 
^  =  t  where  e  is  less  than  1.  Each  of  the  integrals  in  equations 
(B.l)  can  then  be  written  in  the  form 


where  I  is  the  integrand.  The  approximation  (B.4)  is  valid  in  the 

first  integral,  but  not  in  the  second.  However,  provided  X^z  >;>  1, 

_XZ 

the  second  integral  will  be  very  small  because  of  the  e  factor  in 
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the  integrand  and  may  be  neglected.  Thus,  the  condition  for  (B.5)  and 
(B  .6)  to  be  valid  is 

Z  »  T” 

or 


z  »  skin  depth 


Appendix  C 


FOURIER  TRANSFORM  OF  LOOP  CURRENT 


Referring  to  the  coordinate  system  in  Figure  C.l,  the  component 
of  the  loop  current  is  given  by  the  expression 


Jx  =  -I6(z+a)  6(r-R)  sin  0 


(C.l) 


and  therefore  its  Fourier  transform  is  given  by 


J  - - |r  dr  d0  dz  6(z+a)  6(r-R)  sin  0  exp  j  [kAr  cos(0-X)+k  z] 

XK  “  Z 

8jt 


(C .  2) 


where  (r,  0,  z)  and  (k^  X,  kz)  are  polar  coordinates  for  £  and  Jk  and  6 
is  the  Dirac  delta  function.  Write  sin  0  as 


sin  0  =  sin  (0  -  X)  cos  X  +  cos  (0  -  X)  sin  X 


(C.3) 


and  note  from  symmetry  that  the  first  term  will  give  no  contribution  to 
the  integral.  Integrating  the  resulting  equation  yields 


2njR  Ik 

’*k  =  '  sAjl 


eXP  •  (C.4) 


In  similar  fashion 


yR 


2«jR  Ik 

_ x 

8n3k  R 


J1(kj-Ro)  exP  (“Jkza> 


(C.5) 


Equations  (C.4)  and  (C.5)  lead  immediately  to  the  result 


2jkxM 

4  -  -  73—  Ji<k‘V  exp  ('JV> 

8n  k.R 


(C,  6) 
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where 


M  -  jtR  I  i 


(C.7) 


and  i  is  the  unit  vector  in  the  z-direction,  Only  the  case  where 
~z 

k^R  «1  is  of  interest,  and  therefore 


4  = " 


j  kxM 

s*3 


exp  (-j  k  a) 
z 


(  C,  8) 


is  obtained  as  the  final  result. 
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